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A  MULTIVARIABLE  REGRESSION  ALGORITHM 
by 

George  L.  Blaisdell  and  Todd  Carpenter 

INTRODUCTION 

It  is  commonly  the  case  in  data  analysis  that  a  set  of  experimental 
data  is  best  represented  by  an  equation.  The  equation  allows  for  conve¬ 
nient  interpolation  between  measured  data  points  and  provides  a  simple  way 
of  referring  to  the  data.  Additionally,  some  mathematical  manipulation  of 
data  requires  an  equation  form. 

Many  algorithms  exist  for  fitting  a  smooth  curve  to  a  given  i et  of 
data.  It  is  generally  true,  however,  that  these  algorithms  are  limited  to 
a  single  type  of  mathematic  form.  For  instance,  curve  fitting  routines  for 
polynomial,  exponential  and  logarithmic  functions  are  abundant.  These 
curve  types  may  produce  an  adequate  fit  over  a  portion  or  perhaps  all  of 
the  experimental  data.  If  no  one  type  of  mathematical  relationship  ade¬ 
quately  fits  the  data,  one  can  resort  to  using  the  cubic-spline  method  for 
interpolating  data  values,  approximating  the  area  integral  over  a  specific 
interval  or  determining  the  first  derivative  at  a  given  point. 

Since  some,  if  not  most,  experimental  data  do  not  strictly  fit  the 
standard  mathematical  relationships  for  which  curve  fitting  routines  exist, 
in  the  past  they  were  either  roughly  approximated  or  left  as  discrete 
data.  To  overcome  this  difficulty,  we  have  developed  an  algorithm  to  allow 
curve  fitting  with  a  user-defined  mathematical  expression.  This  expression 
may  contain  logarithmic,  exponential,  trigonometric  and  other  mathematical 
forms  singly  or  combined.  It  is  also  able  to  handle  up  to  20  separate 
independent  variables.  Using  the  least  squares  method,  the  algorithm 
determines  the  coefficients  for  each  of  the  user-defined  terms  to 
produce  a  best  fit  curve.  Comparisons  of  the  goodness-of-f it  of  various 
mathematical  expressions  to  the  experimental  data  can  be  done  graphically 


with  plots  or  with  a  fit  parameter  generated  by  the  algorithm. 


A  discussion  of  the  theory  and  concepts  used  to  develop  the  algorithm 
and  details  of  the  computer  coding  are  described  in  this  report.  Several 
application  examples  are  also  illustrated. 


THEORY 

The  least  squares  method  of  curve  fitting  is  used  by  the  algorithm  to 
generate  a  best  fit  curve.  This  method  selects  a  fitted  curve  so  as  to 
minimize  the  sum  of  the  squares  of  the  deviations  of  the  given  data  points 
from  the  curve.  In  the  linear  case  (two  variables — one  independent,  one 
dependent)  the  form  of  curve  to  be  fit  is 

y  =  a  x  +  b  (1) 

where 

a  and  b  =  constants  to  be  solved  for 
x  -  independent  variable 
y  =>  dependent  variable. 

A 

Given  a  set  of  data  points,  denoted  by  y^  (where  i  ranges  from  one  to  n, 
the  number  of  given  data  points),  the  deviation  Dj_  between  the  cal  •.elated 
and  given  values  can  be  expressed 

A 

Di  "  yi  -  yi  (2) 

where 

yt  =  a  xt  +  b 


and  x_  is  the  value  of  the  independent  variable  at  the  gi/eu  data  point 

A 

yi- 

Continuing  with  the  least  squares  method  for  the  linear  case,  and 
using  eq  2,  we  are  given  the  squared  deviations  by 


“  ll'j  ~  (a  xt  +  b>] 

The  sum  of  the  squared  deviations,  S,  is  then  defined  as 


S  -  l  i 

i=l 


(4a) 


(4b) 


2 


A  goodness-of-fit  parameter  can  be  defined  from  this  summation  as  an 
average  deviation,  i.e.,  dividing  S  by  n.  To  minimize  the  sum  S,  the 
partial  derivatives  with  respect  to  each  constant  (a  and  b  for  the  linear 
case)  must  be  set  equal  to  zero.  This  results  in 

n  * 

ff  »  I  ~2  xi  [yi  “  (a  \  +  b)]  =  0  (5a) 

i=l 

n  a 

|f  =  l  -2  [7±  -  (a  +  b)]  =0.  (5b) 

By  rearranging  terms ,  eq  5a  and  5b  become 
n  n  n 

I  y,  x  =  a  l  x2  +  b  l  x  =  0  (6a) 

i=l  i=l  i=l 


n  *  n  n 

£y  =ajx  +b£l  =  0  .  (6b) 

i=l  i=l  1  i=l 


Equations  6a  and  6b  represent  a  series  of  simultaneous  equations  that  can 
be  solved  for  a  and  b.  These  constants  can  then  be  put  into  the  equation 
for  the  general  form  of  the  curve  (eq  1)  to  yield  a  best  fit  curve  for  the 


given  values. 

To  expand  the  least  squares  method  to  include  more  than  just  the 
fit  Co  a  straight  line,  a  general  multivariable,  multiform  expression  is 
defined  as 


g  =•  3l  Tj  +  a2  T2  +  ...  +  aN  Tn 


(7) 


where 

g  **  dependent  variable 
aj  ...  «  constants  to  be  solved  for 

T1  ...  Tjj  =  relational  expressions  of  the  independent 
variables 

N  =  number  of  relational  expressions. 

Using  eq  7  to  determine  the  deviation  to  be  summed  and  minimized  yields 

-  [y£  -  tm  +  a2  t2>1  +  ...  +  aN  y,)] 


N 


■  [y±  -  I 

j=l 


a.  T.  ] 
> 


(8a) 
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where  Tj  is  the  jth  relational  expression  of  the  general  equation 
(eq  7)  evaluated  at  the  ith  set  of  known  values  of  the  independent 
variables.  The  sum  to  be  minimized  can  then  be  defined  as 


N 


S  -  l  [yt  - 1  V,  t]  • 

i=l  1  i=l  J  3,1 


(8b) 


j 

Taking  partial  derivatives  with  respect  to  constants  gives 
Sc  n  -  N 

If-  -  -2  I  Tj,  i  [yt  -  I  ai  Ti  i] 

dak  i=l  j=l  j 


(9) 


where  k  varies  between  1  and  N.  Minimizing  eq  3b  by  setting  eq  9  equal  to 
zero  and  simplifying  gives 


N 


1  ^  ~  'i  3-j  ^4  i  ] 

1  =  1  1  j  =  l  J  J’1 


=  0 


(10) 


where  k  =  1,2,3  ...N.  Equation  9  represents  the  general  form  of  the 
partial  derivative  with  respect  to  one  of  the  constants  in  the  regression 
equation  (eq  7).  To  obtain  the  set  of  simultaneous  equations  for  solving 
the  constants  a^,  eq  10  is  written  specifically  for  k  equal  to  l  to  N. 

This  gives  N  equations  and  N  unknowns  and  can  be  solved  using  the  Gaussian 
elimination  method  with  backward  substitution  and  maximal  column  pivoting. 
Before  applying  this  method,  however,  some  algebraic  simplification  is 
advantageous.  We  can  define  a  constant  C  by 


01) 


Then  using  eq  10,  we  can  further  define  by 

N 

rk,i  Xj,i  = 

j=l  i=l  J=i 


C ,  =  r  a.  r  T,  .  T.  =  'v  a.  A,  . 

k  l  j  l  k,i  j  ,i  l  j  k.j 


(12) 


which  produces  a  series  of  simultaneous  equations. 

Although  laborious,  the  Gaussian  elimination  method  with  maximal 
column  pivoting  is  reliable  and  can  be  easily  programmed,  and  it  can  be 
easily  evaluated  by  a  computer.  Coding  ror  the  least  squares  method  is 
contained  in  program  REGRES  (Appendix  A). 
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PROGRAMMING 

The  multiform,  multivariable  algorithm  was  converted  to  programming 
code  for  an  HP  9845B  minicomputer  (user  manuals  referenced  in  Appendix  B). 
The  algorithm  is  readily  adaptable  to  any  computing  system  with  an 
interpreter;  however,  graphics  and  matrix  manipulation  capabilities  should 
be  available. 

For  convenience,  the  algorithm  has  been  broken  into  five  parts  that 
are  each  separately  programmed  (Appendix  A).  The  first  program,  AUTOST, 
provides  for  data  input  (either  from  keyboard  entry  or  reading  from  a  data 
file  already  generated  by  this  set  of  programs)  and  access  to  the  other 
programs  (REGRES,  PLOTS,  EDIT,  LIST).  Program  REGRES  performs  the  actual 
regression  on  the  user-defined  equation.  PLOTS  provides  a  graphical  output 
for  the  data  points  and  a  curve  or  series  of  curves  of  the  calculated 
equation.  The  EDIT  program  is  a  data  editor  that  allows  additions, 
deletions  or  changes  in  the  input  data.  Program  LIST  provides  a  hard  copy 
listing  of  the  input  data. 

Up  to  20  subfiles  can  i>c  defined  at  the  time  of  data  input.  It  is 
important  that  the  data  be  organized  in  these  subfiles  in  a  logical  manner 
(i.e.  holding  all  /ariables  constant  except  one).  All  operations  contained 
in  the  five  programs  can  be  applied  to  the  complete  data  set  cv  any 
combination  of  subfiles  desired. 

Autostart  program 

AUTOST  is  an  observation  definition  program  that  also  provides  access 
to  the  other  programs  in  the  set.  All  keyboard  input  of  data — including 
subfiles,  subfile  names  and  all  observations — can  be  stored  on  the 
assigned  mass  storage  device  in  a  user-defined  filename  that  does  not 
already  exist.  This  data  file  can  be  called  during  future  program  runs  to 
avoid  repeated  keyboard  entry.  Special  function  keys  K0-K6  are  defined  by 
this  program  by  accessing  the  previously  defined  and  stored  key  file  K1 
(file  K1  should  always  exist  on  the  assigned  mass  storage  device).  Special 
function  key  definitions  are  as  follows: 

K0  -  stops  data  input  into  the  current  subfile  and  allows  data  to 
be  put  into  the  subsequent  subfile  if  desired. 

K1  -  calls  the  regression  equation  program  (REGRES)  and  performs 
the  regression. 
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K2  -  calls  the  PLOTS  program  which  plots  the  most  recently  calcu¬ 
lated  regression  with  the  current  data. 

K3  -  calls  the  EDIT  program  where  the  current  data  may  be  modi¬ 
fied  by  adding,  deleting  or  changing  observations  and  add''  • 
subfiles. 

K4  -  stores  the  current  data  under  a  user-defined  filename 
(provided  that  filename  does  not  already  exist). 

K5  -  calls  the  LIST  program  which  lists  all  the  data  or  certain 
subfiles  cn  the  internal  printer. 

K6  -  stops  the  program. 

Key  KO  is  only  valid  during  keyboard  entry  of  data  and  keys  K]  -  K6  are 
only  active  when  the  prompt  "select  K1-K6"  appears  on  the  screen.  An  error 
message  will  result  if  the  keys  are  pressed  at  any  other  time. 

An  outline  of  the  program  flow  for  AUTOST  is  as  follows: 

1.  The  user  is  asked  if  the  special  function  key  definition  prompts 
should  be  suppressed.  If  yes,  go  to  step  3  ("Y"  or  "YES",  CONT);  if  r.o,  go 
to  step  2. 

2.  The  following  list  is  output  on  the  internal  prLnter: 

KO  =  Stop  data  input  for  current  subfile 

K1  =  Regression 

K2  =  Plots 

K3  ■  Edit 

K4  -  Store 

K5  =  List 

K6  -  Quit. 

3.  The  program  asks  if  stored  data  are  being  used.  If  yes,  go  to 
step  11;  if  no,  go  to  step  4. 

4.  Enter  the  number  of  variables  for  the  data  that  are  to  be  input 
from  the  keyboard. 

5.  Enter  the  name  of  the  subfile  that  the  data  are  to  be  stored  under 
(must  be  no  more  than  five  characters  long). 

6.  Input  the  data  for  all  variables  in  the  current  observation.  If 
this  is  the  first  observation  of  the  subfile,  all  variables  must  be 
included  and  separated  by  commas.  After  data  entry  for  the  first  observa¬ 
tion  is  complete,  an  echo  of  the  data  will  be  printed  on  the  screen. 
Following  entry  of  the  first  observation,  the  user  may  indicate  repeated 
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values  of  any  variable  with  a  "+"  sign  (all  variables  must  still  be  sepa- 
rated  by  commas).  It  is  important  that  all  variables  be  satisfied  for  each 
observation,  whether  it  is  a  repeated  value  or  zero.  Dependent  and  inde¬ 
pendent  variables  are  not  differentiated  at  this  point  but  will  be  identi¬ 
fied  later  in  REGRES. 

7.  If  all  data  entry  for  the  current  subfile  is  complete,  press  KO 
and  go  to  step  8;  if  not,  repeat  step  6. 

8.  The  program  requests  a  new  subfile.  If  more  subfiles  need  to  be 
entered  (answer  "yes"  to  the  prompt),  go  to  step  5;  if  no,  go  to  step  9. 

9.  If  :he  data  are  not  to  be  stored,  answer  "no"  to  the  "data 
stored?"  prompt  and  go  to  step  12.  If  the  data  are  to  be  stored,  answer 
"yes"  and  go  to  step  10. 

10.  The  program  requests  a  filename  (one  that  does  not  already  exist) 
and  then  stores  the  data.  Go  to  step  12. 

11.  Enter  the  filename  under  which  the  data  were  previously  stored  by 
this  program. 

12.  Select  key  K1-K6  to  access  the  other  programs. 

Regression  program 

The  REGRES  program  performs  a  least-squares  regression  on  the  data 
defined  by  AUTOST  with  a  multiform,  multivariable  equation  of  up  to  20 
terms  that  is  supplied  by  the  user.  Each  term  is  input  as  a  character 
string  (up  to  70  characters  long)  and  is  defined  by  the  matrix  notation 
X(A,n),  where  A  indicates  the  observation  number  and  n  denotes  the  number 
of  the  independent  or  dependent  variables.  For  example,  a  regression  on 
the  data  in  Table  1  may  be  desired  using  the  form 

y  =  a  ln(x)  +  bx3  +  cx2z  +  d  (13) 

where 

y  =  dependent  variable 
x  and  z  =  independent  variables 
ln(x)  *  first  regression  term 

3 

x  =  second  regression  term 
x2z  =  third  regression  term 

and  a,b,c  and  d  are  regression  constants  and  the  fourth  regression  term  is  1. 


7 


Under  the  format  required  by  REGRES,  eq  13  would  appear  as 

x(A,2)  =  a  *  LOG  (x(A,l))  +  b  *  x(A,l)  **  3 

+  c  *  x(A,l)  **  2  *  x(A,3)  +  d  (14) 

where  the  variables  x,  y  and  z  are  being  denoted  by  the  numbers  l,  2  and  3, 
respectively,  in  the  matrix  notation. 

To  ease  the  input  of  laborious  equations,  an  option  to  define 
intermediate  constants  for  the  regression  equation  was  made  available. 

Such  constants,  called  user-defined  terms,  are  input  as  a  complete  program 
line  so  it  is  important  that  proper  syntax  is  used.  User-defined  terms  are 
advantageous  when  the  regression  equation  contains  repeatable  sections.  In 
the  example  above,  one  could  simplify  the  equation  by  creating  a 
user-defined  term 

N1  =  x(A,  l) 

and  restructing  the  regression  eq  (eq  13)  as 

x(A,2)  =  a  *  LOG(Nl)  +  b  *  N1  **  3 

+  .  *  N1  **  2  *  x(A,3)  +  d.  (15) 

Clearly,  this  example  does  not  show  much  improvement,  but  user-defined 
terms  can  also  be  used  with  any  programmable  statements. 

The  user  is  restricted  in  the  variable  names  that  can  be  used  to  avoid 
interfering  with  existing  variables.  To  be  safe,  Nl,  N2,  N3,  etc.,  are 
recommended  variable  names. 

The  program  flow  for  REGRES  can  be  summarized  as  follows  (characters 
in  brackets  indicate  correct  responses  for  the  example  above): 

1.  If  this  is  not  the  first  time  REGRES  has  been  run  since  power  was 
turned  on,  go  to  step  9;  otherwise,  go  to  step  2. 

2.  The  program  asks  if  the  user  wants  to  define  any  user-defined 
terras.  If  not,  go  to  step  5  ( N] ;  otherwise,  go  to  3. 

3.  Enter  the  number  of  user-defined  terras. 

4.  Enter  each  user-defined  term  as  a  complete  expression  within 
double  quotes. 

5.  Enter  the  number  of  terms  in  this  regression  (4). 

6.  Enter  the  numerical  index  of  the  dependent  variable  [ 2 J . 

7.  Enter,  one  terra  at  a  time,  the  terras  of  the  regression  equation 
within  double  quotes: 
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["LOG  (X(A,1)"]  (then  press  continue) 

["  X(A, 1)  **  3"] 

["  X(A, 1)  **  2  *  X(A,3)"] 

[1]  or  rx(A,l)  **  0"]. 

8.  If  a  program  or  data  file  named  BUFFER  already  exists  on  the 
assigned  mass  storage  medium,  its  name  must  be  changed  or  it  will  be 
deleted  by  the  REGRES  program.  Go  to  step  16. 

9.  Since  REGRES  has  been  run  before,  the  program  still  has  the  pre¬ 

vious  equation  defined.  It  now  asks  if  the  previous  equation  is  to  be  used 
again.  If  yes,  go  to  step  18;  if  no,  go  to  step  10. 

10.  The  program  asks  if  any  of  the  previous  user-defined  terras  are  to 

be  retained.  If  not  go  to  step  13;  otherwise,  go  to  step  11. 

11.  Enter  the  number  of  user-defined  terms  to  be  retained  and  their 
numerical  indicies. 

12.  Enter  any  additional  user-defined  terms. 

13.  The  program  asks  if  any  of  the  terms  of  the  previous  regression 
equation  are  to  be  retained.  If  not,  go  to  step  15;  if  yes,  go  to  step  14. 

14.  Enter  the  number  of  terms  to  be  retained  and  their  numerical 
indicies. 

15.  Enter  the  terras  necessary  to  complete  the  regression  equation,  as 
in  step  7. 

16.  The  program  now  asks  if  all  the  subfiles  are  to  be  used  in  the 

regression.  If  yes,  go  to  step  18;  if  not,  go  to  step  17. 

17.  Enter  the  number  of  subfiles  and  the  subfile  numbers  wanted  for 

the  regression. 

18.  Following  display  of  the  equation  and  the  solved  constants,  the 
program  asks  if  the  user  wants  a  plot  of  the  data  and  the  regression.  If 
no,  the  program  branches  to  step  12  of  AUTOST;  if  yes,  the  program  branches 
to  program  PLOTS. 

Plots  program 

For  aid  in  determining  the  goodness-of-fit  of  the  regression,  a  plot 
of  the  experimental  and  calculated  data  can  be  produced  with  the  PLOTS 
program.  The  program  will  plot  the  dependent  variable  defined  in  REGRES 
against  any  of  the  other  variables.  Two  plotting  options  are  available  to 
attempt  to  clearly  display  the  relationships  of  several  variables  on  one 
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Table  1.  Example  data  for  multivariable  regression  (values  are 
for  variable  y  [n=2]). 


z(n=3) 


x  (n= 

=1)  1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

58 

62 

65 

69 

73 

76 

80 

83 

37 

91 

2 

90 

104 

119 

133 

147 

162 

176 

191 

205 

219 

3 

157 

190 

222 

255 

287 

319 

352 

394 

417 

449 

A 

277 

334 

392 

449 

507 

565 

622 

680 

737 

795 

5 

463 

553 

643 

733 

823 

913 

1003 

1093 

1133 

1273 

6 

731 

860 

990 

1119 

1249 

1379 

1508 

1639 

1767 

1897 

7 

1096 

1272 

1448 

1625 

1801 

1978 

2154 

2330 

2507 

2683 

8 

1573 

1803 

2034 

2264 

2494 

2725 

2955 

3186 

3416 

3646 

9 

2177 

2469 

2760 

3052 

3343 

3635 

3927 

4218 

4510 

4801 

10 

2924 

3284 

3644 

4004 

4364 

4724 

5084 

5444 

5804 

6164 

graph.  To  illustrate  these,  the  data  from  Table  l  are  used.  AIL  of  the 
x,y  combinations  are  shown  plotted  in  Figure  l.  The  first  plotting  option 
displays  the  dependent  variable  on  the  ordinate  with  one  of  the  independent 
variables  on  the  abscissa  and  a  line  representing  a  constant  value  of  a 
second  Independent  variable  as  calculated  by  the  regression  program  (Fig. 
2).  The  second  type  of  display  plots  the  data  points  with  the  dependent 
variable  on  the  ordinate  and  the  chosen  independent  variable  on  the 
abscissa.  It  also  plots  the  calculated  value  from  the  regression  using  all 
of  the  variables  (Fig.  3). 


x 


Figure  1.  Data  from  Table  1,  x  and  y  only. 
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Figure  2.  Program  PLOTS  output  for  example  data 
of  Table  1  (for  z  =  5) ;  the  ordinate  displays  the 
dependent  variable  (y) ,  the  abscissa  displays  an 
independent  variable  (x)  and  the  plotted  curve 
represents  a  constant  value  of  a  second  independ¬ 
ent  variable  (z  =  5). 


S500 
5850 
5200 
4550 
3900 
>  3250 
2600 
1950 
1300 
650 


u*j  lo  r .  co  cti  o 

K 

Figure  3.  Program  PLOTS  output  for  example  data 
of  Table  1  (for  z  =  5) ;  the  ordinate  displays  the 
dependent  variable  (y) ,  the  abscissa  displays  an 
independent  variable  (x) ;  the  *'s  represent  actual 
data  points  and  the  o's  plot  calculated  points 
(from  the  regression)  using  all  of  the  independent 
variables . 


~  (\J  CO  T 

© 
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The  flow  of  program  PLOTS  follows: 

A, 

1.  Define  the  independent  variable  for  the  abscissa  with  its  index 
number. 

2.  If  all  subfiles  are  to  be  plotted,  go  to  step  4. 

3.  Enter  the  number  of  subfiles  and  the  specific  subfile  numbers 
that  are  to  be  used  in  plotting. 

4.  The  program  displays  the  maximum  and  minimum  values  of  the  vari¬ 
able  chosen  for  the  plot  axes.  The  user  is  asked  to  enter  the  limits  for 
the  axes'  scaling. 

5.  Enter  the  label  for  the  abscissa  (independent  variable)  and  the 
ordinate  (dependent  variable). 

6.  To  plot,  holding  all  variables  constant  except  one  (the 
independent  variable  chosen  in  step  1),  go  to  step  11. 

7.  The  program  plots  each  data  point  (*)  and  the  calculated  points 
from  the  regression  program  (0)  for  each  data  point  using  all  variables. 

8.  The  program  asks  if  a  hard  copy  is  desired,  if  not  it  branches  to 
step  10. 

9.  A  hard  copy  is  printed  with  the  regression  equation  and  a 
goodness-of-fit  parameter  (a  fit  parameter  of  0  is  a  perfect  fit).  If  a 
different  plot  is  wanted,  go  to  step  1. 

10.  If  another  regression  is  desired,  the  program  branches  to  step  1 
of  program  REGRES.  If  no  additional  regressions  are  desired,  the  program 
branches  to  step  12  of  AUTOST. 

11.  The  program  requests  entry  of  a  constant  value  for  all  variables 
except  the  one  identified  in  step  1. 

12.  The  plot  is  displayed.  If  another  curve  is  desired,  go  to  step 
13;  if  not,  go  to  step  8. 

13.  The  program  requests  the  number  of  variables  to  be  changed,  which 
ones  (identified  by  their  numerical  index)  and  their  new  values.  Go  to 
step  12. 

Edit  and  list  programs 

Program  EDIT  allows  the  user  to  change  the  values  of  observations,  add 
and  delete  observations,  and  add  complete  subfiles  to  the  existing  files. 
All  editing  features  include  either  subfile  or  absolute  addressing. 

Absolute  addressing  numbers  all  observations  sequentially  from  the  first 
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observation  of  the  first  subfile  to  the  last  observation  of  the  last 
subfile.  Subfile  addressing  renumbers  all  observations  in  each  subfile 
such  that  the  first  observation  of  each  subfile  is  referenced  as  the  first 
observation.  The  EDIT  program  offers  menu-type  edit  feature  selection. 

The  delete  feature  first  lists  the  values  of  all  variables  of  the 
observation  referenced  and  offers  an  escape  by  asking  if  it  should  be 
deleted. 

The  LIST  program  provides  a  hard  copy  listing  of  some  or  all  sub¬ 
files.  It  prints  the  subfile  name  at  the  top  of  each  subfile  and  refer¬ 
ences  all  observation  numbers  using  subfile  addressing. 

ALGORITHM  APPLICATION 

In  the  study  of  snow  mechanics,  one  of  the  common  techniques  used  to 
classify  snow  compressibility  is  the  plate-sinkage  test.  A  circular  or 
rectangular  plate  driven  by  a  hydraulic  or  manual  ram  compacts  a  volume  of 
snow,  changing  its  original  density  to  critical  density  (approximately 
3.12x10  lb/' ft  ).  At  the  critical  density,  the  volume  of  snow  is  reduced 

very  little  or  not  at  all.  Increased  pressure  from  the  plate  results  in 
constant-volume  flow  or  movement  of  the  snow  away  from  the  plate. 

Measurements  taken  during  a  plate-sinkage  test  include  the  force  on 
the  plate  and  an  indication  of  the  amount  of  sinkage  or  vertical  compaction 
experienced  between  the  beginning  of  the  test  (zero  force)  and  each  force 
reading.  From  the  force-sinkage  data,  the  compaction  energy  can  be  deter¬ 
mined  by  integration  between  the  zero-sinkage  and  critical  sinkage  (sinkage 
at  which  critical  density  is  reached)  limits.  This  integration  can  be 
performed  numerically  to  yield  acceptable  results.  Integration  of  the 
curve  or  comparison  of  the  nature  of  the  compaction  relationship  between 
various  snows,  however,  is  accomplished  best  with  a  mathematical  represen¬ 
tation.  Using  the  multivariable  algorithm,  we  can  determine  the 
mathematical  relationship  between  force  and  sinkage. 

Given  the  force  versus  sinkage  data  in  Table  2,  a  best-fit  mathemati¬ 
cal  expression  is  desired.  Since  we  know  that  the  curve  is  somewhat 
exponential  in  nature,  the  first  regression  equation  form  used  in  REGRES  is 

F  =  a  eZ  (16) 

which  gives  aj  =  0.01347  and  from  PLOTS  shows  a  rather  poor  fit  (Fig.  4). 
The  second  equation  type  attempted  is  a  4th  order  polynomial  and  it  results 
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Table  2.  Typical  force  versus  sinkage  data 
for  a  plate-sinkage  test  in  snow. 


Force  (lb) 


Sinkage  (in.) 


0 

10 

15 

15 

20 

25 

35 

45 

50 

70 

90 

115 

150 

170 

200 

240 

275 

320 

355 

400 

460 

535 

610 

830 


0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 
4.0 


4.5 
5.0 

5.5 

6.0 

6.5 
7.0 

7.5 

8.0 

8.5 

9.0 

9.5 

10.0 

10.5 

11.0 

11.5 


in  the  coefficients  and  fit  shown  in  Figure  5.  This  is  a  much  better 
approximation  to  the  data,  however,  near  the  upper  and  lower  ends  of  the 
data  the  polynomial  shows  an  increasingly  poor  fit.  (This  is  a  common 
behavior  for  polynomial  curve  fitting.)  Recognizing  that  the  exponential 
example  showed  a  horizontal  curve  near  the  low  sinkage  values  and  became 
steeper  at  the  high  sinkage  values,  we  next  attempt  a  combination  of 
exponential  and  polynomial  forms.  The  equation  used  by  REGRES  for  this 
trial  is 

«  n  k  2 

F  =  ai  +  &2Z  +  a3z  +  a4Z  +  asz  +  a6e  .  (17) 
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SIMKRGE  (inches) 
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»■-.  1315148723155 
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<  — 14. 1377573413 
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N2"EXP<X<A,2>> 

Goodn«ii-»f  fu-tt  43.825 


Figure  6.  Fourth-order  polynomial  plus  exponential 
curve  fit  for  plate-siukage  data. 


Performing  the  regression  with  eq  17  and  the  data  in  Table  2,  we  obtain  an 
excellent  fit  (Fig.  6).  Another  form  of  combining  the  polynomial  and 
exponential  fits  is  shown  in  Figure  7  where 

F  =  ez  (a!  +  a2z  +  a3z2  +  a4z3  +  a5z4)  +  a5  .  (18) 

This  also  results  in  a  good  curve  fit  and  comparison  of  the  goodness-of-f it 
parameter  is  necessary  to  determine  that  eq  18  more  closely  approximates 
the  data. 

Typical  heat  transfer  problems  demonstrate  the  multivariable  power  of 
this  regression  algorithm.  A  case  in  point  involves  the  placement  of  flat 
heat  flux  sensors  on  a  curved  surface,  such  as  a  pipe.  Heat  flux  sensors 
(HFS)  are  thin  wafers  of  a  material  with  known  thermal  conductivity. 
Thermocouples  are  attached  to  each  flat  side  of  the  wafer.  The  output  from 
the  HFS  is  a  voltage  that  is  proportional  to  the  heat  flux  through  the 
sensor.  The  flat,  inflexible  HFS  is  attached  to  the  pipe  by  an  epoxy.  The 
epoxy,  in  addition  to  securing  the  sensor  to  the  pipe,  fills  in  the  void 
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)<|"X<A,2> 

H2-EXP<X<A,2>> 

Goodn* s»-of  fU-»»  20.337 

Figure  7.  Fourth-order  polynomial  times  exponential 
curve  fit  for  plate-sinkage  data. 

space  between  the  flat  HFS  and  the  curved  pipe.  This  attachment  alters  the 
factory  calibration  for  thermal  conductivity,  which  is  for  application  on  a 
flat  surface.  Using  experimental  data  from  an  insulated  pipe  section  with 
a  HFS  attached,  we  can  determined  a  new  calibration  for  thermal  conduc¬ 
tivity  by  relating  the  heat  flux  out  of  the  pipe  to  the  sensor  /oltage 
output.  This  calibration  relation  can  easily  be  represented  by  a  smooth 
curve.  Different  insulation  thicknesses  generate  a  series  nf  curves. 
Additionally,  with  different  pipe  insulation  thermal  conductivities,  pages 
of  these  series  of  curves  are  generated.  Likewise,  with  different  glues 
for  attaching  the  sensor,  volumes  of  pages  of  series  of  curves  are 
required.  It  soon  becomes  apparent  that  a  mathematical  relation  is 
necessary  to  calibrate  the  heat  flux  out  of  the  pipe  with  the  HFS  reading, 
the  type  of  glue,  the  type  and  thickness  of  the  insulation,  and  the 
diameter  of  the  pipe. 

The  first  step  in  fitting  an  equation  to  these  data  is  to  break  the 
data  into  subfiles  where  only  one  variable  is  changing.  This  isolates  the 
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effect  each  variable  has  on  the  heat  flux  out  of  the  pipe  and  the  heat  flux 
through  the  sensor.  Starting  with  the  basic  equation  for  heat  flux  out  of 
the  pipe 


♦ 


P 


(19) 


where 

%  =  insulation  thermal  conductivity 

<J>p  =  heat  flux  out  of  pipe 

•Pd  =  heat  flux  through  sensor 

fp  =  radius  of  the  pipe 

tj[  =  thickness  of  insulation, 

we  determine  the  effect  of  each  variable  on  the  heat  flux.  All  the  terms 


HERT  FLUX  THRU  SENSOR  i  BTU  Hr-  ft  2  F> 
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Figure  8.  Multivariable  curve  fit  for  heat  flux 
through  an  insulated  pipe  with  a  flat  heat  flux 
sensor. 
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of  these  equations  can  then  be  included  in  a  final  regression  on  all  the 
subfiles.  The  final  form  of  the  best  fit  equation  xs 


♦ 


P 


aK,  <j>  dK 

—  —  -7  +  b  K  <j>  +  cK  +  — | 
r  +  t.  irs  i  <J> 


1/3 


In 


i  K  <f> 

+  f  <b  K  +  g  r  t,  +  h  K  r  t.  + - =— -  +  j  (20) 

yss&pi  spi  r  t.  J 

pi 


where  Ke  is  the  epoxy  thermal  conductivity  and  Ks  is  the  sensor  thermal 
conductivity.  The  constants  a-j  and  a  plot  of  the  calculated  and  raw  data 
are  shown  in  Figure  ".  This  figure  appears  cluttered  because  all  the 
subfiles  of  data  were  plotted.  This  can  be  avoided  by  either  expanding  the 
scale  to  look  at  discrete  sections  of  the  data  or  plotting  one  subfile  at  a 
time. 


CONCLUSION 

The  multivariable,  multiform  regression  algorithm  presented  can  be 
seen  to  be  a  useful  alternative  to  single  form  regressions.  The  algorithm 
allows  data  that  do  not  fit  one  of  the  standard  equation  forms  (polynomial, 
exponential,  etc.)  to  be  mathematically  characterized  by  a  best-fit  curve. 
This  allows  the  data  to  be  easily  integrated  and  differentiated  and 
provides  an  accurate  interpolating  equation.  An  equation  representation  of 
the  data  also  provides  a  method  of  comparing  the  variable  relationships 
between  separate  sets  of  data. 

Specific  application  of  the  algorithm  is  shown  for  snow  compaction 
with  a  plate-sinkage  device.  The  results  can  be  used  for  predicting 
vehicle  motion  resistance  when  the  vehicle  is  operating  in  snow.  The 
algorithm  is  also  applied  to  the  problem  of  heat  flow  around  a  pipe. 
Measured  heat  flow  data  can  be  fit  with  a  regression  equation  and  a 
comparison  of  actual  and  analytical  solutions  completed. 
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APPENDIX  A:  ALGORITHM  CODING. 


10  REM  — >  PROGRAM  AUTOST 

20  !  STARTER  PROGRAM  FOR  MULTI-VARIABLE  REGRESSION. 

30  PRINT  PAGE,  TP-.B<  13)  J  CHR*  <  129);  "USER  INSTRUCTIONS  FOR  MULT  I -VARIABLE  FEGFE33 
ION: " ; CHRSC 12S) 

40  PRINT  LINC2); "This  program  is  the  starter  program  for  a  set  of  five  progra 
ms" 

50  PRINT  "that  can  perform  a  series  regression  of  a  user-aefined  equation  mt 
h  " 

SO  PRINT  "up  to  500  observations  and  20  variables  per  observat i on. " 

70  PRINT  LINCl);"The  data  is  input  from  the  terminal  and  can  later  be  stored 

for  future  use." 

80  PRINT  "All  data  is  broken  into  subfiles  that  can  be  any  length  so  long  as 
the  500  " 

90  PRINT  "observation  maximum  is  not  exceeded.  Up  to  20  subfiles  can  be  spec 
i f i ed, " 

100  PRINT  "and  any  regression  or  plotting  function  can  be  performed  with  »n»  c 
onibinat  ion" 

110  PRINT  "of  these  subfiles." 

120  PRINT  LINCl);"Any  variable  can  be  assigned  the  dependent  variable  by  dssig 
n&t  i  ng" 

130  PRINT  "the  index  number  of  the  variable  desired,  therefore  all  variables  a 
re  treated" 

140  PRINT  "the  same.  Plots  can  be  made  with  any  variable  as  the  independent  v 
ar  i  ab I e  " 

150  PRINT  “except  the  dependent  variable.  The  dependent  variable  can  not  appe 
ar  in  the" 

160  PRINT  "regression  equation." 

170  DISP  "PRESS  'CONT'  TO  CONTINUE." 

180  PAUSE 

190  PRINT  PAGE, TABC30);CHR*C132>; "PROGRAM  0PTI0NS";CHRS- 128) 

200  PRINT  LlNC2>;“T!ie  following  programs  and  options  ar-/  available!" 

210  PRINT  LINC 1 ) , TABc5) { “  1  >  Keyboard  data  input  with  subfile  definition,  or  re 
ad  data  from" 

220  PRINT  TAB<8>J"a  previously  stored  data  file." 

230  PRINT  LIN< 1 ) , TAB<5> ; "2)  Regression  of  user  defined  series  equation  such  as 
Y=a*X''3+b*l  n<X)+c  “ 

240  PRINT  TAB<8> J "where  a,b,and  c  are  the  constants  that  the  program  solves." 

250  PRINT  L INC  1 ) , TAB<5> 5 "3>  Plot  the  regression  equation  with  a  data  scatter  p 

lot  to  show  correlation" 

260  PRINT  LINC 1 ) , TABC5)  } "4 )  Data  editing  including  observation  changes,  additi 
ons,  and  deletions  with" 

270  PRINT  TFiB<8) ;  "subf  1 1  e  additions  and  data  file  merging." 

280  PRINT  LINC 1 ) , TABC5) J "5>  Store  the  current  data  in  a  user-defined  filename. 

41 

290  PRINT  LINC 1 ) , TABC5) ; “6>  List  all  the  data  or  only  selected  subfiles.' 

3O0  DISP  "PRESS  'CONT'  TO  CONTINUE." 

313  PAUSE 

320  PRINT  PAGE,  TABC22)  JCHR^C  132)  J  "DEFINITIONS  OF  USER-DEFINED  t-  E  ■  3“ ;  LHP4  -  1  23  » 

330  PRINT  LINC2),TABc5>; "The  user-defined  keys  K0-K6  are  used  to  quiet  1'  acces 

s  the  various" 

340  PRINT  “program  options.  K0  is  used  during  data  input  to  suspend  input 
the  current'1 

350  PRINT  “subfile.  It  is  used  only  after  all  data  has  been  recorded  into  the 
current  " 

360  PRINT  "subfile.  The  option  of  adding  a  subfile  is  then  made  available,  a* 
ter  which" 

370  PRINT  "all  data  is  put  into  the  new  subfile." 
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333  PRINT  LIN< 1 } , ThB<5> ; "The  keys  K1-K6  access  the  o*  her  progr  am  options,  and 
should  be  used  " 

330  PRINT  "only  when  the  prompt  SELECT  K1-K6  appears.  The  plot  option  shoul 
d  not  be" 

400  PRINT  "used  before  a  regression  has  been  performed  since  it  does  not  ha"s 


410  PRINT  "equation  to  perform. " 

420  PRINT  LIN< 1 >, TAB<5) ; " The  option  of  getting  a  hard  copy  reminder  of  what  th 
ese  keys  do" 

430  PRINT  "is  now  available,  you  can  suppress  the  prompts  if  so  desired." 

440  OPTION  BASE  1 

450  DIM  A*<20>tl23,B*<20>C12] 

460  COM  X<500, 20) , V$ <20>  C 160] , Nsubf <20; C 103 , Abt 20, 21 ) , INTEGER  Dep, Subf • 20  ■ . Nsu 
b, Nobs , Nv, Nt , M$<20) C 150  3 , Term 


470 

480 

490 

500 

510 

520 

530 

540 

550 

551 
5?0 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 

770 

771 
730 
790 
800 
810 
820 
830 
340 
350 
360 
370 
330 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
AT 


INPUT  "SUPRESS  KEY  PROMPTS?" , S* 

G0SUB  Check 
IF  Check=l  THEN  470 
IF  UPC*<S*Ci,  13>="Y"  THEN  590 
PRINTER  IS  0 

PRINT  LIN<2> , "K0  =  STOP  DATA  INPUT  FOR  CURRENT  SUBFILE." 
PRINT  "K1  =  REGRESSION.  " 

PRINT  "K2  =  PLOT." 

PRINT  "K3  =  EDIT.  " 

PRINT  " K4  =  STORE.  " 

PRINT  "K5  =  LIST  " 

PRINT  "K6  a  QUIT" 

PRINTER  IS  16 
GCLEAR 

EXIT  GRAPHICS 
LORD  KEY  "Kl:T15" 

PRINT  PAGE 

INPUT  "ARE  YOU  USING  STORED  DATA?",S$ 

GOSUB  Check 
IF  Check* 1  THEN  640 
IF  UPC*<S*C1, 13>="N"  THEN  770 
ON  ERROR  GOTO  Badfile 
INPUT  "FILE  NAME'’"  ,  F  i  1  nm$ 

ASSIGN  til  TO  Fi  1  nm$ 

OFF  ERROR 
PRINT  PAGE 

READ  # 1 ; Nv, Nobs, Nsub 

REDIM  X<Nobs, Nv) , Subf <Nsub+ 1 ) , Nsubs <Nsub> 

READ  #l;Nsub$<*>,  Subf  <  *  ) , X<  *  > 

GOTO  1360 

INPUT  "HUMBER  OF  VARIABLES?  <\  =  20V,  Nu 
REDIM  X<500, Nv) 

IMAGE  #,2X,4A,2D 
Hsum=0 

IMAGE  #,1X,7A 
IMAGE  #, 1X.5D.4D 
Nobs*Subf <  1  )=0 
Nsub=0 

Nsub=Nsub+ 1 

DISP  "NAME  OF  SUBFILE  #";Nsub;“  = 

INPUT  "  <<=  5  CHARACTERS > " , Nsub* ( Ns ub > 

PRINT  CHR*<27)&"m" 

PRINT  PAGE, TAB<30); "  SUBFILE  -  " ; Nsub*  < Nsub  * 

PRINT  "CBS#" 

FOR  B=1  TO  Nv 

PRINT  USING  780;"VAR#“;B 
NEXT  B 
PRINT 

PRINT  CHR*  < 27 )  & "  1  11 
REDIM  A*<Nv>,B*<Nv> 

Bcount=0 

Nobs=Nobs+ 1 
Bcount=Bcount+l 

DISP  "INPUT  ALL  VARIABLES  FOR  OBSERVATION  # " ; Be ount ; " 

II  « 

> 


REPE 
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1000  MAT  INPUT  A* 

1910  IF  BcountOl  THEN  1050 

1020  FOR  C=1  TO  Nv 

1030  B*<C)=A*<C> 

1040  NEXT  C 

1050  PRINT 

1060  PRINT  LIN<2> ; Bcount 

1070  FOR  C=1  TO  Nv 

1080  IF  A$<C ><>"+"  THEN  1140 

1090  GOTO  1170 

1100  X(Nobs, C)=VAL(B$<C>  ) 

1110  A*<C)=B*<C> 

1120  PRINT  USING  810; XOlubs, C> 

1130  GOTO  1170 

1140  IF  A*<C>="+"  THEN  1170 

1150  X<Nobs,C>=VAL<A*<C>> 

1160  PRINT  USING  810; X<Nobs, C> 

1170  B*<C>=A*<C> 

1180  NEXT  C 

1190  GOTO  970 

1200  Swbfi! INPUT  "ANOTHER  SUBFILE?" ,Sf 

1210  GOSUB  Check 

1220  IF  Check=l  THEN  SubFi 

1230  Nobs=Nobs- 1 

1240  SubF<Nsub+l )=Nobs 

1250  PRINT  CHR*<27)$<"m"  ,  PAGE 

1260  IF  L'PC$(S9C 1 , 1 3  >■" Y"  THEN  840 

1270  INPUT  "STORE  DATfi^'.S* 

1280  IF  UPC*<S*C1, U)»"N"  THEN  1360 

1290  Phyrec  =  <Nohs*20*12''256+Nsub*?2',256+2>*3-'2 

1300  REDIM  A*<Nv) , B#<Nv) , Nsubf <Nsub> , Subf <N3ub+ 1 > , XOIobs , Nv> 

1310  INPUT  "FILE  NAME?" , F l 1 nm* 

1320  CREATE  F l 1 nm* , Phyrec 
1330  ASSIGN  #1  TO  Fi lnm* 

1340  PRINT  H 1 ; Nv, Nobs , Nsub 

1350  PRINT  HI ; Nsub$<* J , Subf <*) , X< * ) 

1360  Restart:  Flaggg*5© 

1361  OH  KEY  HI  GOTO  1450 
1370  ON  KEY  H2  GOTO  1470 
1380  ON  KEY  H3  GOTO  1490 
1390  ON  KEY  H4  GOTO  1290 
1400  ON  KEY  H5  GOTO  1510 
1410  ON  KEY  H6  GOTO  1440 

1420  DISP  "SELECT  K1-K6" 

1421  Flaggg=Flaggg+l 
1430  GOTO  1430 

1440  END 

1450  DISP  "GOING  TO  REGRESSION  PROGRAM." 

1460  LOAD  "REGRES: T15" 

1470  IF  FlagggM  THEN  GOTO  1479 

1471  DISP  "PLOTS  CAN  ONLY  BE  ACCESSED  AFTER  REGRESSION  HAS  BEEN  RUN 

1472  WAIT  6000 

1473  DISP  “SELECT  K1-K6" 

1474  GOTO  1430 

1479  DISP  “GOING  TO  PLOTTING  PROGRAM." 

1480  LOAD  “PLOTS! T15“ 

1490  DISP  “GOING  TO  EDIT  PROGRAM." 

1500  LOAD  “EDIT:T15" 

1510  DISP  "GOING  TO  LISTING  PROGRAM." 

1520  LOAD  "LIST: T15" 

1530  Check:  Check=0 
1540  PRINT  PAGE 

1550  IF  (UPC*<S*t 1 , 1]>  =  "Y",>  OR  <UPCf<S*Cl, i]>=“N";  THEN  1590 
1560  Check=l 
1570  BEEP 

1580  PRINT  PAGE;"***  IMPROPER  RESPONSE  -  TRY  AGAIN  ***" 

1590  RETURN 

1600  !  CHECKS  FILENAME 

1610  Badf 1 le:  BEEP 
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1620  IF  ERRN056  THEN  1650 

1630  PRINT  PAGE;"-***  FILE  NAME  DOES  NOT  EXIST  **+" 
1640  GOTO  690 
1650  DISP  ERRMi 
1660  GOTO  1440 


10  REM  — >  PROGRAM  REGRES 

20  OPTION  BASE  1 

30  COM  X<500, 20) , V$<  20) C  1 60  3 ,  Nsubi<20>  C103,Ab<.20,21),  INTEGER  Dep ,  3ubf  ■  20  ' ,  Ni  u 
b, Nobs , Mo, Nt , M*<20) C 1 503 , Term 
32  REDIM  X<Nobs, No) 

40  DIM  Dp#C1603,St$tl6O3,StliC1603,S)<20>,Keep<2O),Ni<2G>U5O3,Test  fC  1503.  Prf 
C  1503 

50  INTEGER  User<20> 

60  MAT  Ab=2ER 

70  Slct=l 

80  Kterm=0 

90  DISP 

100  Last=N=Nt 

110  IF  Vi < 1 ) = " "  THEN  370 

120  Last  =N=Nt 

140  INPUT  "USE  SAME  EQUAT I  ON?" , Si 

150  GOSUB  Check 

160  IF  Check=l  THEN  140 

170  IF  UPCi<SiC 1 , 1 3  >  =  "Y"  THEN  2340 

180  INPUT  "DO  YOU  WANT  TO  KEEP  ANY  OF  THE  TERMS  OF  THE  EGUAT I  ON-" , Si 

190  GOSUB  Check 

20O  IF  Check=l  THEN  180 

210  IF  UPC«<S$C1, 13)a“N"  THEN  670 

220  INPUT  "HOW  MANY  TERNS  DO  YOU  WANT  TO  KEEP^",Knt 
230  IF  (Knt  >0)  AND  <Knt<=20>  AND  <FRACT (Knt  >=0)  THEN  270 

240  BEEP 

250  PRINT  PAGE; "***  IMPROPER  NUMBER  OF  TERMS  ***" 

260  GOTO  220 

270  PRINT  PAGE 

280  FOR  A=1  TO  Knt 

290  INPUT  "TERM  #?",Keep(A) 

300  IF  (Keep(A) >0)  AND  (Keep(A)< =20)  AND  > FRACT(Keep( A  * )=0 >  THEN  ?40 

310  BEEP 

320  PRINT  PAGE;"***  IMPROPER  TERM  NUMBER  ***" 

330  GOTO  290 

340  PRINT  PAGE 

350  NEXT  A 

360  GOTO  670 

370  PRINT  PA6E,TAB*.25<;CHRi(129>;“N0TES  ON  REGRESSION  PROGRAM" ;  CHRi  'IIS' 

380  PRINT  LIN( 2) , TABt5 * ; “The  regression  program  takes  an  equation  defined  b  t 

he  user  such  as!" 

390  PRINT  LIN<  1>;  “  Y=a*X~3+b*X*  1  n(2>+c  " 

400  PRINT  LIN<l);"and  soloes  for  a,b,and  c  based  on  the  best  fit  to  the  curren 
t  data.  “ 

410  PRINT  L INt 1 ) , TAB( 5) ; " The  relational  expressions  are  the  input  terms  of  *  he 
regression  equation." 

420  PRINT  "Once  you  haoe  input  the  regression  equation,  you  can  perform  the  " 

430  PRINT  "regression  on  any  combination  of  subfiles  to  find  the  relation  in  o 

ne  " 

440  PRINT  “particular  subfile.  Later,  eoen  after  a  plot  has  been  performed,  ■ 
ou  c  an  “ 

450  PRINT  "repeat  the  same  regression  without  the  need  of  retyping  the  reg-ess 
i on  equat i on" 
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460  PRINT  "or  the  user  defined  terms.  So  you  can  perform  the  regression  on  on 
ly  a  feu" 

470  PRINT  "subfiles  the  first  time,  and  later  perform  the  regression  on  »he  "h 
ole  data  set . " 

480  PRINT  "Also,  you  can  pick  and  choose  which  regression  terms  and  which  user- 
defined  " 

490  PRINT  "terms  you  wish  to  keep,  so  remember  the  order  in  which  you  input  th 
e  terms  to  " 

5O0  PRINT  "capitalize  on  this  feature." 

510  DISP  "PRESS  'CONT'  TO  CONTINUE." 

520  PAUSE 

530  PRINT  PAGE, TAB<20);CHRT<129);" HON  TO  INPUT  USER  DEFINED  TERMS" J CHRT ■  1 28  1 
540  PRINT  LIN<2) , TAB<5) ; "User  defined  terms  are  relational  expressions  that 
complement  the  " 

550  PRINT  "regression  equation  by  defining  constants  that  can  be  t  •  ?d  in  ti-.;  e 
quat ion." 

560  PRINT  "These  terms  are  independent  equations  t h a ,  must  have  proper  syrt  a  . 


and  any" 

570  PRINT  "variable  names  can  be  used  except  A,B,  and  any  variable  in  common. 
User  defined  " 

580  PRINT  "terms  are  placed,  in  order,  before  the  regression  equation,  and  are 
eval uated" 

590  PRINT  "for  every  point  in  the  data  set.  It  is  important  that  the  s,,nt a  o 
f  any  user" 

600  PRINT  "input  terms  is  correct.  If,  for  any  reason,  the  program  gets  lost" 
610  PRINT  "<i.e.  the  screen  is  blank,  and  the  tape  drive  is  not  working1  for  a 
n  extended" 

620  PRINT  "length  of  time,  stop  the  program,  and  start  over.  The  problem  was 
due  to  " 

630  PRINT  "improper  syntax  of  the  input  equations." 

640  PRINT  LINC 1 ), TAB<5) ; "User  defined  terms  can  also  be  used  with  any  programa 
ble  function" 


650  PRINT  "such  as  PRINT  or  RAD  to  display  any  intermediate  results  or  tempora 
ri ly  set  the" 

660  PRINT  "computer  into  a  desired  computational  mode." 

670  IF  NS < 1 ) a " "  THEN  880 

680  INPUT  "DO  YOU  NANT  TO  KEEP  ANY  OF  THE  USER  DEFINED  TERMS*" , ST 

690  GOSUB  Check 

700  IF  Check=>l  THEN  680 

710  IF  UPC*<StCl, 13)="N"  THEN  880 

720  INPUT  "HON  MANY  USER  DEFINED  TERMS  DO  YOU  NANT  TO  KEEP'1" , Kt erm 
730  IF  <Kterm>0)  AND  <Kterm<20)  AND  <FRACTfKterm >=0)  THEN  770 
740  BEEP 

750  PRINT  PAGE; IMPROPER  NUMBER  OF  USER  DEFINED  TERMS  ***" 

760  GOTO  720 

770  PRINT  PAGE 

780  FOR  A= 1  TO  Kterm 

790  INPUT  "KEEP  USER  DEFINED  #",UsertA) 

800  IF  (Usrr  <A)  >0)  AND  < User  <  AX  =20 >  AND  «  FRACTtUser  <.  A)  <  =  0  >  THEN  340 

810  BEEP 

820  PRINT  PAGE;"***  IMPROPER  USER  DEFINED  TERM  NUMBER  ***" 

830  GOTO  790 

340  PRINT  PAGE 

350  N«<User<A))=NT<User>.A)) 

360  NEXT  A 

370  GOTO  930 

330  INPUT  "DO  YOU  NANT  ANY  USER  DEFINED  TERMS'*"  ,  ST 

390  GOSUB  Check 

900  IF  Check=l  THEN  880 

910  Term=0 

920  IF  UPC*<S*C1, 1]>="N"  THEN  1140 

930  INPUT  "HOW  MANY  USER  DEFINED  TERMS'*" ,  Term 

940  IF  <Term>0)  AND  tTerm<=20>  AND  <FRACT >  Term)=0>  THEN  980 

950  PRINT  PAGE;"***  IMPROPER  NUMBER  OF  USER  DEFINED  TERMS  **♦" 

960  BEEP 

970  GOTO  930 

980  PRINT  PAGE 

990  FOR  A=1  TO  Term 
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1000  FI ag=0 

10ie  FOR  B=1  TO  Kterm 

1020  IF  User (B)=A  THEN  Flag=l 

1030  NEXT  B 

1040  IF  FlagOl  THEN  1070 

1050  N*(A)=M*(A) 

1060  GOTO  1090 

107O  DISP  "INPUT  USER  DEFINED  TERM  #";A;"  UN  QUOTES;  150  character  ma  lmuw 

II  • 

1080  INPUT  " " , N$( A) 

1090  Test$=NS(A) 

1100  GOSUB  Parenth 

1110  IF  Ptest=l  THEN  1070 

1120  M*(A)=N$(A) 

1130  NEXT  A 

1140  PRINT  PAGE ,  TAB<  13,1 ;  CHRf  <  129 )  J  "PROCEDURE  FOR  THE  INPUT  OF  THE  REGRESSION  ECl 
UAT ION" ; CHR$  ( 12S) 

1150  PRINT  LIN<4) ; CHRit 132) ; "FOR  THE  EQUATION  ' X ( A , 3 ,  =  I *X < A , 2 >  2  +  J-X' A, 1  '  - 
K' “ J  CHR$( 128) 

1160  PRINT  L I N < 1 ) ; "  X  ( A ,  3  )  =  THE  DEPENDENT  VAPIABLE.  IN  THIS  CASE  THE  THIRD  V 

ARI ABLE. " 

1170  PRINT  LINO); "  I 4X< A, 2>A2  =  THE  FIRST  TERM  OF  THE  EQUATION  FOR  WHICH  I 
IS  TO  BE  SOLVED." 

1180  PRINT  "  (Input  'X(A,2)'2'  as  the  first  term  within  double  quotes.  »" 

1190  PRINT  LINO);"  J-*X(A,1>  =  THE  SECOND  TERM  FOR  WHICH  J  IS  TO  BE  SOLVED." 

1200  PRINT  LINO);"  K  =  THE  THIRD  TERM  WHICH  IS  A  CONSTANT  WITHOUT  A  VARIABLE 

II 

1210  PRINT  "  (Input  O'  to  represent  this  term.)" 

1220  PRINT  LIN( 1 ); "Regression  terms  are  a  maximum  of  70  characters  in  length." 
1230  Start: INPUT  "HOW  MANY  TERMS  ARE  THERE  IN  THE  REGRESS  I  ON Nt 
1240  IF  (Nt  >0)  AND  (Nt<»20>  AND  (FRACT ( Nt ) =0 )  THEN  1280 
1250  BEEP 

1260  PRINT  PAGE;"***  IMPROPER  NUMBER  OF  TERMS  ***" 

1270  GOTO  Start 
1280  PRINT  PAGE 
1290  REDIM  Ab ( Nt , Nt  + 1 ) 

1300  PRINT 
1310  Last  =N=Nt 

1330  INPUT  "INPUT  THE  INDEX  OF  THE  DEPENDENT  VARIABLE,  ve.g.  1  the  for  1st  wari 
bl e)“ , Dep 

1340  IF  (Dep>0)  AND  (Dep<.=Nv)  AND  (FRACT ( Dep>=0)  THEN  1380 
1350  BEEP 

1360  PRINT  PAGE;"***  IMPROPER  VARIABLE  NUMBER  ***" 

1376  GOTO  1330 
1380  PRINT  PRGE 

1390  V*(Ht  +  l  )  =  "X<'A,  "S<VAL*(Dep)4<")" 

140O  FOR  A= 1  TO  Nt 

1410  FI ag=0 

1420  FOR  B= 1  TO  Knt 

1430  IF  A=Keep<  B  >  THEN  Flag=l 

1440  NEXT  B 

1450  IF  Flag=l  THEN  1590 

1460  DISP  "INPUT  TERM  #";A;“  (e.g.  X»'A,U  within  quotes*"; 

1470  INPUT  11 "  ,  Vt  (  A ) 

1480  Test*=V*(A) 

1490  GOSUB  Parenth 

1500  IF  Ptest=l  THEN  1460 

1510  St=POS(V$(A),V*(Nt+l>) 

1520  IF  St=0  THEN  1590 

1530  PRINT  USING  “+, K" ; Vf i A > ; “  IS  NOT  VALID  SINCE  THE  DEPENDENT  VARIABLE  I: 
";V*(Nt+l) 

1540  BEEP 

1550  DISP  "PRESS  CONT  TO  CONTINUE" 

1560  PAUSE 

1570  PRINT  USING  "  +  ,K";  " 

II 

1580  GOTO  1460 

1590  NEXT  A 
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1600  PRINT 

1610  ON  ERROR  GOTO  1630 
1620  GOTO  1680 
1630  IF  ERRN054  THEN  1660 
1640  PURGE  "BUFFER: T15" 

1650  GOTO  1680 
1660  PRINT  ERRM* 

1670  PAUSE 

1680  CREATE  "BUFFER: T15" , <Nt-KNt+Nu+Term+3>+5>*S0/256+l 
1690  ASSIGN  "BUFFER: T15"  TO  #1 
1700  PRINT  #l;"50  Builder:!" 

1710  PRINT  #  1 ; "51  FOR  B=1  TO  Fx" 

1720  Ln=50 

1730  PRINT  H  1 ;  VAL$<51  +Ln)!!<"  FOR  A=Subf  <  S 1  <  B  >  )  + 1  TO  SubfCSl <B>+1 V 

1740  Ln=Ln+l 

1750  FOR  A=1  TO  Term 

1760  PRINT  Hl;VAL*<51+ln>8<"  "«<N$(A) 

1770  Ln=Ln+l 

1780  Pol=POS<N*<A>, "X<A, "> 

1790  IF  Pol=0  THEN  Poend 

1800  Po2=Pol +POS<N*< A)  CPol  ] ,  "  ) "  ) 

1810  N$<A)CPol,  Po2-l  ]  =  "S1  < "&N$ < A) [Pol +4 ,  Po2-2H"  )  " 

1820  GOTO  1780 

1830  Poend!  NEXT  A 
1840  FOR  R=1  TO  Nt 
1850  FOR  B=1  TO  Nt+1 

1860  St*=VAL$<51+Ln)ij<"  Ab<  "8,VALi<  A>«<"  ,  "S<VALi<B>8."^  =  "S<"Pb-  "S<VRLT*A< 

ALf<B>&")  +  "$iV$<A)&"*"&V$(B) 


1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

I960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2O80 

2090 

21O0 

2110 

2120 

2130 

2140 

2150 


LnsLn+1 
PRINT  Hi;St$ 

NEXT  B 
NEXT  A 

PRINT  H 1 J  VAL* <51 +ln)&"  NEXT  A" 

PRINT  HI ;  VAl* (52+LnH"  NEXT  B" 

PRINT  Hl;VAl*(53+LnH"  GOTO  5500  " 
ln<*Ln+3 
FOR  D*1  TO  Nu 
FOR  B=1  TO  Nt 

St*="  Fx=Fx+"  W<BH"-*Ab<"&VAl*<BH",  "kVALJOit  +  l  >8,"  >" 
St  U=St* 

Strt=POS<St 1$, "X<A,  "&VAL$<D)  > 

IF  Strt=0  THEN  2040 
Nd=POS<St lJCStrt  ],")"> 

St  UCStrt  ,Strt+Nd-n  =  "Var " 

GOTO  1990 
FOR  C= 1  TO  Nv 

Str=POS(St  1 1,  " X C  A ,  "iVALKC)  > 

IF  Str=0  THEN  2100 
En=POS<Stl*CStr],">"> 

St  l*CStr,  Str+En-ll  =  "Sl  <  “e.VAL*<  C  >&“  >" 

GOTO  2050 
NEXT  C 

IF  BOl  THEN  2240 

PRINT  HI ;  VAL*<51+LnH"  A'&VALJCDH":  Fx=0" 

Ln=Ln+l 

FOR  Ed=l  TO  Term 
Pr*=N*<Ed> 


2160  Pol=POS<Pr*,  "SlC'fcVAL^D)) 

2170  IF  Pol=0  THEN  2210 

2188  Po2=Pol +POS<Pr$ [Pol 

2190  Pr$CPol,Po2]="Var" 

2280  GOTO  2160 

2210  PRINT  H1JVAL*<51+Ln>*."  "8cPr* 

2220  Ln=Ln+ 1 

2230  NEXT  Ed 

2240  PRINT  Hl;VAL*<51+Lnn"  "8<Stl* 

2250  Ln=Ln+l 

2260  NEXT  B 

2270  PRINT  Hl;VAL$<51+Ln>«."  GOTO  5500" 
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2280  Ln=Ln+l 

2290  St  li=St* 

2300  NEXT  D 

2310  PRINT  #l;"5508  SUBEND" 

2320  ASSIGN  #1  TO  * 

2330  OFF  ERROR 

2340  LINK  "BUFFER:T15", 3450, 2350 
2350  PRINT  PAGE 

2360  INPUT  "DO  YOU  WANT  TO  DO  THE  REGRESSION  OH  ALL  SUBFILES'5" , S$ 

2370  GOSUB  Check 

2380  IF  Check=l  THEN  2360 

2390  IF  UPC$<S*C1, 1])="Y"  THEN  2560 

2400  INPUT  "HOW  MANY  SUBFILES  DO  YOU  WANT  INCLUDED  IN  THE  REGRESS  I  ON?" , Fx 
2410  IF  <Fx>0)  AND  <Fx<=Nsub)  AND  <FRACT<Fx>=0>  THEN  2460 

2420  DISP  "***  WARNING  -  NUMBER  OF  SUBFILES  EXCEEDED  +  **" 

2430  BEEP 
2440  WAIT  2000 

2450  GOTO  2400 

2460  FOR  A=1  TO  Fx 

2470  DISP  "INPUT  SUBFILE  #";A;"  FOR  THE  REGRESSION."; 

2480  INPUT  " " , S 1 <  A  ) 

2490  IF  <S1  <A>  >0)  AND  <SKAX=Nsub>  AND  <FRACT<S1  <A>  >  =  0>  THEN  2530 
2500  BEEP 

2510  PRINT  PAGE;"***  IMPROPER  SUBFILE  NUMBER  ***" 

2520  GOTO  2470 

2530  PRINT  PAGE 

2540  NEXT  A 

2541  PRINT  PAGE, LIN(8),TAB<25): rHRf < 129) ; "REGRESSION  IN  PROGRESS" : CHRIU 28) 
2550  GOTO  2610 

2560  FOR  ft= 1  TO  Nsub 
2570  S1(A)=A 

2580  NEXT  A 

2590  PRINT  PAGE,  L  INO)  ,  TAB05)  ;  CHRS <  129)  ; "REGRESSION  IN  PROGRESS" ; CHRX'  1 2S  1 
2600  Fx=Nsub 

2S10  CALL  Bui ld<S1ct,Var,Fx,Sl <*),N$<*)) 

2620  !  SIMULTANEOUS  EQUATION  PROGRAM  OF  Nth  ORDER. 

2630  FOR  Mat  =  l  TO  N-l 

2640  First=Mat 

2650  Large=First 

2660  FOR  Test=First+l  TO  Last 

2670  IF  Ab< Test , F irst > >Ab(Large , F i rst )  -THEN  Large=Test 

2680  NEXT  Test 

2690  IF  Large=F i rst  THEN  GOTO  2750 

2700  FOR  Switch=First  TO  Last+1 

2710  Save=Ab<Fi rst , Swi tch) 

2720  Ab (First, Suite h)=Ab<Large, Switch) 

2730  Ab<Large, Swi tch)=Saoe 

2740  NEXT  Switch 

2750  FOR  Diag=First+l  TO  Last 

2760  Rat i o=Ab<  D i ag , F i rst )/Ab<F irst, First) 

2770  FOR  2ero=First  TO  Last+1 

2780  Ab<Di  ag, Zero)=Ab<  Di  ag, Zero^-Rat i o*Ab  (First,  Zero  > 

2790  NEXT  Zero 

2800  NEXT  Di ag 

2810  NEXT  Mat 

2820  FOR  Sol ue=Last  TO  1  STEP  -1 

2830  Ab<Sol  ve ,  L  ast  + 1  i=Ab  t  So  1  we ,  L  ast  +  1  XAb<$ol  ve,  Sol  ve  > 

2840  IF  So  1 oe=  1  THEN  GOTO  2890 

2850  FOR  Sol v=$o 1 ve- 1  TO  1  STEP  -1 

2860  Ab<Solo,Soloe)=Ab<Soloe,Last+l)*Ab<Solo,Solue) 

2870  Ab<Sol  u,  Last +  1  )=Ab<Sol  o,  Las t  + 1 ) -AbOo  1  o,  Sol  oe  > 

2880  NEXT  Sol v 

2890  NEXT  Solve 
2900  PRINT  PAGE 
2910  PRINT  LINO) 

2920  PRINT  USING  " t* ,  k“  ;  "Fx=a*"  ;  V*<  1 ) 

2930  FOR  8=2  TO  Nt 

2940  PRINT  USING  " # , K" ;"+"; CHRS <96+R );"+"; V$ < A > 
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2950  NEXT  ft 

2960  PRINT 

2970  FOR  fl=l  TO  Nt 

2980  PRINT  TflB<3) ;  CHRS (96+ft ) ; "  =  " ; flb<fl, Nt  +  1  > 

2990  NEXT  ft 

3000  FOR  fl=l  TO  Term 

3010  PRINT  LIN<l);M$<fl) 

3020  NEXT  ft 

3021  BEEP 

3022  WAIT  1000 

3023  BEEP 

3030  INPUT  "DO  YOU  WANT  ft  PLOT?",S* 

3040  GOSUB  Check 

3050  IF  Check=l  THEN  3030 

3060  IF  UPC*<S*C1, n)="Y"  THEN  LOAD  "PLOTS ! T 1 5" 

3070  PRINT  PAGE 

3080  DISP  "RETURNING  TO  MAIN  PROGRAM" 

3090  LOAD  "AUT0ST:T15", Restart 
3100  ! 

3110  !  PARENTHESIS  CHECKING  ROUTINE 

3120  Parenth:  Ptest=0 

3130  Left  count=Right  count =Poss=Po=0 
3140  Po=POS<Test$, " < 

3150  IF  Po=0  THEN  Leftpas 
3160  Test$CPo, Po3="*“ 

3170  Left  count=Left  count+1 
3180  GOTO- 3140 

3190  Leftpas:  Poss=POS<Test t , "  > " ) 

3200  IF  Poss=0  THEN  Rightpas 
3210  Test*[Poss,Poss]="*“ 

3220  Right  count=»Right  count  +  1 
3230  GOTO  Leftpas 

3240  Rightpas:  IF  Left_count=Ri ght_count  THEN  RETURN 
3250  Ptest=l 

3260  IF  Left  count>Right  count  THEN  DISP  Left  count-Pight  count ; "EXCESS  -ErT  PA 
RENTHESIS" 

3270  IF  Left  count<Right  count  THEN  DISP  Right  count-Left  count ; "EXCESS  RIGHT  P 
ARENTHESIS" 

3280  BEEP 
3290  WAIT  3000 
3300  RETURN 
3310  ' 

3320  !  CHECKS  YES  AND  NO  RESPONSES 

3330  Check:  Check=0 
3340  PRINT  PAGE 

3350  IF  <UPCI<S*tt,ll>«"Y">  OR  <UPC*<SfC 1, 1J>»"N“>  THEN  RETURN 
3360  Checks  1 
3370  BEEP 

3380  PRINT  PAGE;"***  IMPPOPER  RESPONSE  -  TRY  AGAIN." 

3390  RETURN 
3400  ! 

3410  SUB  Bui ld<Slct,Var,Fx,Sl <*>,N$<*>> 

3420  OPTION  BASE  1 

3430  COM  X<*>, V*<*),Nsub$<*>, Ab<*>, INTEGER  Dep, Subf <* » , Nsub, Nobs , Nu, Nt 
3440  ON  SI ct  GOTO  Buil der , A t , A2, A3, A4 , A5, A6, A7, A8, A9, A1 0, A 1 1 , A1 2, ftl 3, A1 4 . A 1 5 . A1 
6, A17, A18, A19.A20 
3450  Bui lder: ! 


10  REM  — >  PROGRAM  LIST 

20  OPTION  BASE  1 

30  COM  X<500, 20)  ,  V$<20>  C 1603 , Nsubif  20  )  C 103 , ftbi 20, 2 1  ■ ,  INTEGER  Dep.Subf'20  .tisu 
b, Nobs, No, Nt , M$<20>  C 1583 , Term 
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32  REDIM  X<Nobs,Nv) 

49  DIM  Sub<20> 

50  INPUT  "DO  YOU  WANT  ALL  SUBFILES  LISTED?", S$ 

so  if  <upc$<s*ci,  n)="Y">  or  <upc$<s$ci,  m  =  "N">  then  si 
70  BEEP 

SO  PRINT  PAGE; "444  IMPROPER  RESPONSE  -  TRY  AGAIN  tt." 

90  GOTO  50 

91  INPUT  "DO  YOU  WANT  A  HARD  COPY  <H>  OR  A  CRT  DISPLAY  'OV.Tf 

92  IF  < UPC#< T$  [  1 , 1 ] ) ="H " >  OR  <UPC* < T* C 1 , 1 ] >  =  "C " ;  THEN  96 

93  BEEP 

94  PRINT  PAGE;"  444  IMPROPER  RESPONCE,  PLEASE  TRY  AGAIN  4-4" 

95  GOTO  91 

96  IF  UPC$<T$C1, 1 ] ) = "C"  THEN  101 

97  PRINTER  IS  0 

101  PRINT  PAGE 

110  IF  UPC$<S$C1, 1])="Y"  THEN  230 

120  INPUT  "HOW  MANY  SUBFILES  DO  YOU  WANT  L  I STED'’ "  ,  Nss 
130  IF  <Nss>0)  AND  <Nss<=l'sub>  AND  <FRACT<Nss>=0>  THEN  170 
140  BEEP 

ISO  PRINT  PAGE; "444  IMPROPER  NUMBER  OF  SUBFILES  *'4" 

160  GOTO  120 

170  PRINT  PAGE 

180  FOR  A=1  TO  Nss 

190  DISP  "INPUT  SUBFILE  #";A; 

200  INPUT  " " , Sub< A) 

210  IF  <Sub<A>>0>  AND  < Sub< AX  =Nsub >  AND  <FRACT<,Sub'.A>  ;=0  »  THEN  250 

220  BEEP 

230  PRINT  PAGE; "444  IMPROPER  SUBFILE  NUMBER  «4»" 

240  GOTO  190 

250  PRINT  PAGE 

260  NEXT  A 

270  GOTO  320 

230  FOR  A= 1  TO  Nsub 

290  Sub< A)=A 

300  NEXT  A 

310  Nss=Nsub 

340  FOR  R*1  TO  Nss 

341  PRINT  LIN<3> 

350  PRINT  L  IN<  1  >  ;  "SUBF  ILEl  ",NsubF<.A> 

360  PRINT  LIN( 1 ) 

370  Jj=0 

380  FOR  I i=l  TO  Nv 

390  Jj=Jj+l 

400  IF  Jj<=5  THEN  GOTO  430 

410  Jj=i 

420  PRINT  USING  460 

430  PRINT  USING  450;  “VAR#"S,VAL*<  1 1  > 

440  NEXT  li 

450  IMAGE  #,7X,6A 

460  IMAGE  / 

470  PRINT  USING  460 

480  PRINT  "  OBS  #" 

490  Obs=0 

500  FOR  B=Subf <$ub<A) )+l  TO  Subf <Sub< A)+l > 

510  Obs=Obs+l 

520  PRINT  LINC1) 

530  PRINT  Obs 

540  Jj=0 

550  FOR  I i=l  TO  Nv 

560  Jj=Jj+1 

570  IF  Jj<=5  THEN  GOTO  60O 

580  Jj=l 

590  PRINT  USING  460 

600  PRINT  USING  680;X<B,li> 

610  NEXT  Ii 

620  NEXT  B 

630  NEXT  A 

631  PRINT  L I N< 2 ) 
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640  PRINTER  IS  16 

650  PAUSE 

660  DISP  "RETURNING  TO  MAIN  PROGRAM" 

670  LOAD  "AUT0ST:T15", Restart 

680  IMAGE  #,SX,5D.4D 

690  IMAGE  10X 


10  REM  — >  PROGRAM  PLOTS 

20  OPTION  EASE  1 

30  COM  X<500, 20) , Vt<20) C 160] , Nsubf<20> C 103 , Ab<20, 21 >, INTEGER  Dep, Subf <. 20 > .  Hsu 
b, Nobs, Nv, Nt , M*<20) C  150] , Term 
32  REDIM  X<Nobs,Nv) 

40  DIM  Gxmxx<20>  ,  Gxmnx'  20) ,  Gymxx<28>  ,  Gymnx<  20)  ,  S 1  t  20.) ,  Eqf  C500] ,  XI  a*  C 1 60  ] ,  VI  as 
C 160] 

60  INTEGER  Nf  i le<20) 

70  PRINT  PAGE, TAB <25 > ; CHR* < 1 29 ) ; "NOTES  ON  PLOTT I NG: " ; CHRS <  1 28) 

30  PRINT  L IN<2) , TAB< 5) ;  “Thi  s  program  will  plot  t he  data  and  the  regression  fu 

notion  using  any" 

90  PRINT  "variable  as  the  independent  variable  except  the  dependent  variable. 


100  PRINT  LIN< 1 > , TAB<5) J "The  function  can  be  plotted  using  a  line  and  setting 
al 1  var i abl es" 


110  PRINT  "to  a  constant  except  the  independent  and  dependent  variables.  Mul» 
i pi e  1 i nes" 

120  PRINT  "can  be  drawn  in  this  manner  with  the  option  of  changing  any  of  the 
constants. " 


130  PRINT  "Another  option  is  to  plot  correspond i ng  data  points  that  niser's  al 
1  variables" 

140  PRINT  "from  each  data  point  into  the  equation." 

150  PRINT  LIN< l ) , TAB<5) ! "Any  combination  of  subfiles  can  be  plotted." 

160  LINK  "BUFFER:T15", 2340, 170 

170  GCLEAR 

180  EXIT  GRAPHICS 

190  DEG 

200  INPUT  "INDEPENDENT  VAR  I ABLE?" , Ind 
210  IF  IndODep  THEN  240 

220  PRINT  PAGE,LIN<7); "THAT  IS  THE  DEPENDENT  VARIABLE." 

221  GOTO  200 

230  GOTO  200 

240  PRINT  PAGE 


250  INPUT  "PLOT  THE  COMPLETE  DATA  SET>",S* 

260  GOSUB  Check 

270  IF  Check=l  THEN  250 

280  IF  UPC*<$*C1,  n>  =  ,,N“  THEN  380 

290  Sbsts=Nsub 

300  MAT  SEARCH  X<*, Ind) , MIN; Xmn 
310  MAT  SEARCH  X< * , Ind) , MAX; Xmx 
320  MAT  SEARCH  X< * , Dep) , M IN; Ymn 
330  MAT  SEARCH  X<*, Dep) , MAXJ  Ymx 
340  FOR  A= 1  TO  Nsub 
350  Nf l 1 e<A)=A 

360  NEXT  A 

370  GOTO  650 

380  INPUT  "HOW  MANY  SUB  F I LES?" . Sbs t s 

390  IF  <Sbsts>0)  AND  <Sbsts<=Nsub J  AND  <FRACT<Sbsts)=0)  THEN  430 
400  PRINT  PAGE!'1***  IMPROPER  NUMBER  OF  SUBFILES  ***" 

410  BEEP 

420  GOTO  380 

430  PRINT  PAGE 

440  REDIM  GxmnxtSbst  s ) , Gxmxx<Sbst  s  > , GymnxCSbst s ) , GymxxtSbst  s  > 

450  FOR  A=1  TO  Sbsts 
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•*60  INPUT  "NUMBER  OF  SUBF I LE? “  ,  Nf  1  1  e ( A> 

470  IF  <Nf i 1 e<RP  >0)  AND  < Nf l 1 e < A > < =Nsub>  AND  (FRACT <Nf i 1 e< A  >  >=0>  THEN 

480  BEEP 

480  PRINT  PAGE;  "-<•*■*  IMPROPER  SUBFILE  NUMBER 

500  GOTO  460 

510  PRINT  PAGE 

520  Gxnmx<  A)  =Gxmxx<  A>=X<Subf  <Nf  i  I  «<  A)  )  + 1 ,  Ind) 

530  Gymnx< A>=Gymxx< A)=X<Subf <Nft 1 e< A> >+l , Dep> 

540  FOR  B=Subf<Nfi  1e<'A>)  +  l  TO  Subf <Nf l 1 e< A>+ 1 > 

550  IF  X<B,  I nd ) >Gxmxx < A >  THEN  Gxrnxx<A>- XCB,  Ind> 

560  IF  X<B,  IndXGxmnx<A>  THEN  GxrnnxC A>  =X( B ,  I nd > 

570  IF  X< B, Dep) >Gymxx< A)  THEN  Gytoxx<A)=X<B, Dep> 

580  IF  X(B,  Dep)<Gyninx<  A>  THEN  Gymnx<A>=X<B, Dep> 

590  NEXT  B 

600  NEXT  A 

610  MAT  SEARCH  Gxmxx, MAX; Xmx 

620  MAT  SEARCH  Gxmnx, MIN; Xmn 

630  MAT  SEARCH  Gymxx, MAX; Ymx 

640  MAT  SEARCH  Gymnx, MIN; Ymn 

650  PRINT  PAGE;"X-MIN  =  " ; DR0UND<Xmn,3> ; "  XHAX  =  " ; DROUND<Xmx , 3 > 

660  PRINT  L I N < 2 ) ; " Y-M I N  =  “ ; DROUND< Ymn , 3 > ; "  YMAX  =  " ; DROUND*  Ym  ,3< 

670  INPUT  "INPUT  GRAPH  LIMITS  FOR  X-MIN  AND  X-MAX. " , Gxmri, Gxfftx 

680  INPUT  "INPUT  GRAPH  LIMITS  FOR  Y-M IN  AND  Y-MAX . " , Gy  dm , Gvmx 

690  PRINT  PAGE 

700  LIMIT  0, 184.47,0, 149.3 

710  LOCATE  15,  120, 15,95 

720  CSI2E  3.38 

730  SCALE  Gxtnn,  Gxmx,  Gymn,  Gymx 

740  FRAME 

750  DISP  "HANG  ON  A  MINUTE" 

760  Xt ica<Gxmx-Gxmn)/10 
770  Ytic  =  <Gy(iix-Gytnr>)/'10 
780  AXES  Xt )  c  ,  Yt  i  c  ,  Gxrrm,  Gymn 
790  LORG  8 
800  LDIR  90 

810  FOR  A=Gxmn  TO  Gxmx  STEP  Xtic 

820  MOVE  A, Gymn 

830  LABEL  DR0UND<A,3> 

840  NEXT  A 

850  INPUT  "X-AXIS  LABEL", XI a« 

860  DISP  “HANG  ON  A  MINUTE" 

870  LDIR  0 

880  FOR  A-Gymn  TO  Gymx  STEP  Ytic 

890  MOVE  Gxmn.A 

900  LABEL  DROUND<  A , 3> 

910  NEXT  A 

920  LORG  6 

930  MOVE  <Gxmx+Gxmo)^2 , Gymn-<Gyrox-Gymn > '8 

940  LABEL  XI  a* 

950  INPUT  "Y-AXIS  LABEl ", Y 1  at 

960  MOVE  Gxmn-<Gxmx-Gxrrin)/8.  1 ,  <Gymx+Gymn>/2 

970  LORG  4 

980  LDIR  90 

990  LABEL  Ylal 

1000  LORG  5 

1010  INPUT  "DO  YOU  WANT  AN  EQUATION  LlNEtl;  OR  CORRESPONDING  DATA  Pi'iINT*. 
1020  IF  <S=0>  OR  <S=1)  THEN  1060 
1030  BEEP 

1040  PRINT  PAGE;"***  IMPROPER  RESPONSE  -  TRY  AGAIN  **■*" 

1050  GOTO  1010 
1060  PRINT  PAGE 
1070  Slct=Ind+l 
1030  IF  S=0  THEN  1430 
1090  Pais=0 

1100  IF  PaisO0  THEN  1120 
1110  GOTO  1370 

1120  INPUT  "DO  YOU  WANT  THE  SAME  PARAMETERS'’" ,  it 

1130  GOSUB  Check 

1140  IF  Check=l  THEN  1120 
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1158  IF  UPC*<S*tl,  1]>  =  ,,Y"  THEN  1430 

1160  INPUT  "HOW  MANY  PARAMETERS  DO  YOU  WANT  TO  CHANGE?", S 
1170  IF  <S>0>  AND  <S<=Nw?  AND  < FRACT<$ > =0 >  THEN  1210 
1180  PRINT  PAGE;"***  IMPROPER  NUMBER  OF  PARAMETERS  ***" 

1190  BEEP 
1200  GOTO  1160 
1210  PRINT  PAGE 
1220  FOR  A=1  TO  S 

1230  INPUT  "WHAT  VARIABLE  NUMBER?" , Vnn 

1240  IF  <Vnn>0>  AND  <.Vnn<=Nu>  AND  <FRACT<Vnn>=0>  THEN  1280 
1250  BEEP 

1260  PRINT  PAGE;"***  IMPROPER  VARIABLE  NUMBER  ***" 

1270  GOTO  1230 

1280  PRINT  PAGE 

1290  IF  <VnnOD*p>  AND  (VnnOInd)  THEN  1320 

1300  PRINT  PAGE, LIN(7>; "THAT  IS  EITHER  THE  INDEPENDENT  OR  THE  DEPENDENT  VAP 
ABLE.  " 

1310  GOTO  1230 

1320  PRINT  PAGE 

1330  DISP  "WHAT  VALUE  DO  YOU  WANT  FOR  VARIABLE  #";Vnn; 

1340  INPUT  " " , S 1 (Vnn) 

1350  NEXT  A 
1360  GOTO  1420 
1370  FOR  A= 1  TO  Nu 

1380  IF  <AaInd>  OR  <A=Dep>  THEN  1410 

1390  DISP  "WHAT  VALUE  DO  YOU  WANT  FOR  VARIABLE  t*";A; 

1400  INPUT  " " , S 1 < A ) 

1410  NEXT  A 

1420  Pass  = 1 

1430  GRAPHICS 

1440  Vari  anee=Nbcount.=0 

1450  FOR  A-l  TO  Sbsta 

1460  FOR  B=Subf<m'i  1e<A>>-M  TO  Subf CNf 1  I *< A>  +  1 ) 

1470  Nbcount=Nbcount+l 

1480  MOVE  X<B, Ind),X<B, Dep> 

1490  LABEL  "*" 

1500  Var=X<B,Ind> 

1510  FOR  C=»l  TO  Nu 

1520  IF  <C»Dep>  OR  <C»Ir,d)  THEN  1540 

1530  SI <C)=X<B,C) 

1540  NEXT  C 

1550  CALL  Bui  1  d  <  S 1 c  t , Var, Fx , SI <  *  )  > 

1560  Var i ance=Var i ance+<Fx-X<B, Dep>  >A2 

1570  IF  SO0  THEN  1600 

1530  MOVE  Var , Fx 

1590  LABEL  "0" 

1600  NEXT  B 

1610  NEXT  A 

1620  IF  S<  >0  THEN  1710 

1630  LORG  2 

1640  MOVE  Gxmn+Xt i c^3, Gymx-Yt i c^3 
1650  CSI2E  2.7 
1660  LDIR  0 

1670  LABEL  "0  =  CALCULATED" 

1680  MOVE  Gxmn+Xt  l  c  '3,  Gymx"3/'4*Yt  i  c 
1690  LABEL  “*  =  DATA" 

1700  GOTO  1840 

1710  FOR  R=Gxmn  TO  Gxmx  STEP  <Gxmx-Gxmn)/50 
1720  Var=A 

1730  CALL  Bui  ld<Slct, Var.Fx.Sl <*)) 

1740  IF  A=Gxmn  THEN  MOVE  Var,Fx 

1750  PLOT  Var, Fx 

1760  NEXT  A 

1770  PAUSE 

1780  EXIT  GRAPHICS 

1790  INPUT  “ANOTHER  LINEP",S* 

1800  GOSUB  Cheek 

1810  IF  Check=l  THEN  1790 

1320  IF  UPC*<S*C1, n?=“Y"  THEN  1120 
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1830  GOTO  1850 
1840  PAUSE 

1850  INPUT  "DO  YOU  WANT  A  HARD  COPY?",S* 

I860  GOSUB  Check 

1870  IF  Check=l  THEN  1850 

1880  IF  UPC$<S*C1,  13)="N"  TKS.N  2050 

1890  PRINTER  IS  0 

1900  PRINT  PAGE 

1910  DUMP  GRAPHICS 

1920  Eq$="Fx=a*"S,V$0) 

1930  FOR  A=2  TO  Nt 

1940  Eq*  =  Eq*V  +  "&CHR$<9t  < A) 

1950  NEXT  A 
1960  PRINT  Eq* 

1970  FOR  A= 1  TO  Nt 

1980  PRINT  CHR*<A  +  96)$."  =  "8<VAL*'"2b<A,Nt  +  n) 

1990  NEXT  A 

2000  FOR  A= 1  TO  Term 

2010  PRINT  LINO);M*<A> 

2020  NEXT  A 

2030  PRINT  LINO);  "Goodness-of  fit  is  " ;  PROUNDc  V  ar  i  ance'  Nbcount ,  -  3 ) 

2040  PRINTER  IS  16 

2050  INPUT  "ANOTHER  PLOT?" , Si 

2060  GOSUB  Cheek 

2070  IF  Check=l  THEN  2050 

2080  IF  UPC*<S*C1, 1I)="Y"  THEN  170 

2090  PRINT  PAGE 

2100  GCLEAR 

2110  EXIT  GRAPHICS 

2120  INPUT  "ANOTHER  REGRESS  I  ON?" , S* 

2130  GOSUB  Check 

2140  IF  Check=l  THEN  2120 

2150  IF  UPC*<S*C1,  m  =  "Y"  THEN  DISP  "GOING  TO  REGRESSION  PROGRAM." 

2160  IF  UPC*<S*C1,  1 3  >  = "  Y  "  THEN  LOAD  "RE'.'RES!  T  1 5" 

2170  DISP  "RETURNING  TO  MAIN  PROGRAM" 

2180  LOAD  "AUTOSTJTIS", Restart 
2190  ! 

2200  !  CHECKS  YES  AND  NO  RESPONSES 

2210  Check:  Checks© 

2220  PRINT  PAGE 

2230  IF  <UPC*<S*C 1 , 1  ]  >  =  "Y" )  OR  <UPCI <$* C 1 , 1 3 >  =  "N" >  THEN  2270 
2240  BEEP 
2250  Check=l 

2260  PRINT  PAGE;"***  IMPROPER  RESPONSE  -  TRY  AGAIN  ***" 

2270  RETURN 
2280  1 
2290  ■ 

23O0  SUB  Bui ld<Slet ,Var,Fx,S) (*)) 

2310  OPTION  BASE  1 

2320  COM  X<*),V$i.*),Nsub$<  +  ),AbO>,  INTEGER  Dep,  Subf  <.  *  > ,  Nsub,  Nobs ,  Nu,  Nt  ,  Mf 
rni 

2330  ON  Slct  GOTO  Bu i 1 der , A l , A2, A3, A4 , A5, A6, A7, A8 . A9, A1 0, A1 1 , A1 2, A1 3. A1 4 , A I 

6,A17,A18,A19,A20 

2340  Bui  1 der:  ! 


10  REM  — /  PROGRAM  EDIT 

20  OPTION  BASE  1 

30  COM  XC500, 20) , Vf i 20  1 C 1601 , Nsubf i 20^ C 101 , Abf 20. 21  i , INTEGER  Dep.Subf'20 
b, Nobs , Nw, Nt , M$s  20) 1 1 501 , Term 

31  REDIM  X(500,Nu) 

40  Menu:  PRINT  PAGE, LIN< 3); "SELECT  AN  EDITOR  FUNCTION:" 

50  PRINT  LIN<1)5"  1)  CHANGE  AN  OBSERVATION." 

60  PRINT  LINO);"  2)  CHANGE  A  SET  OF  OBSERVATIONS." 
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70 
80 
80 
100 
1  10 
120 
130 
140 
150 
180 
170 
180 
180 
2O0 
210 
220 
230 
240 
250 
280 
270 
280 
280 
300 
310 
320 
330 
340 
350 
380 
370 
330 
380 
400 
410 
420 
430 


PRINT  LIN<1); 
PRINT  LIN<1>; 
PRINT  LIN<1); 
PRINT  LIN<1>; 
PRINT  LIH< 1 ) ; 
"  "  ,  Ef 
AND 


3)  ADD  A  SUBFILE. " 

4)  ADD  AN  OBSERVATION." 

5)  DELETE  AN  OBSERVATION." 
8>  MERGE  A  DATA  FILE." 

7) 


<Ef<8>  AND  <FRACT<Ef>=0>  THEN  ISO 


QUIT  THE  EDITOR." 

INPUT 
IF  < Ef >0 > 

BEEP 

PRINT  PAGE;"***  IMPROPER  RESPONSE  ♦**" 

WAIT  3000 
GOTO  Menu 
PRINT  PAGE 

ON  EF  GOTO  20O,  Chgset ,  Add, Addo , De 1 ete , Merge , Ou i K 


INPUT  "DO  YOU  WANT  SUBFILE  ADDRESSING  (O'*  OR  ABSOLUTE  ADDRESSING  •  1  • '",3 

IF  <S=0>  OR  <S=1)  THEN  280 

BEEP 

PRINT  PAGE;"***  IMPROPER  RESPONSE  ***" 

GOTO  210 

PRINT  PAGE 

IF  S=1  THEN  Absol 

INPUT  "WHICH  SUBFILE  #?",Sn 

IF  <Sn>0)  AND  <Sn<=Hsub>  AND  <FRACT<.Sn>=Q »  THEN  330 
BEEP 

PRINT  PAGE;"***  SUBFILE  NUMBER  HOT  DEFINED  ***" 

GOTO  280 
PRINT  PAGE 

DISP  "WHICH  OBSERVATION  IN  SUBFILE  " ;  lUubf <Sn- ; 

INPUT  Subn 

IF  <Subn>0)  AND  <$ubn<=$ubf <Sn+l »-Subfv$n>>  AND  1 FRACT < Subr  '=0  >  THEN  400 
BEEP 

PRINT  PAGE;"***  SUBFILE  OBSERVATION  ADDRESS  OUT  OF  RANGE  ►**" 

GOTO  340 
PRINT  PAGE 

INPUT  "WHICH  VARIABLE  tt" ,  Vn 

IF  (Vn>0)  AND  <Vn<=Nw  AND  < FRACT ( Vn ) =0 )  THEN  480 
BEEP 


440  PRINT  PAGE!"***  IMPROPER  VARIABLE  NUMBER  ♦**" 

450  GOTO  410 
480  PRINT  PAGE 

470  DISP  "OLD  VALUE  ■  “  ;  X<SubF<Sn>  *Subn,  Vn)  ,*  "  NEW  VALUE'"! 

480  INPUT  , A 

480  PRINT  USING  510;  "OBS#! " , Subn, "  IN  SUBFILE!  ",  NsubT<Sn; ,  "  OLD  VALUE!  ",  Su 

bF<Sn>+Subn,Vn), "  NEW  VALUE!", A 

500  X<Subf<Sn)+Subn, Vn)=A 

510  IMAGE  5A,4D, 13A.5A, 12A.5D.5D, 12A,5D. 5D 

520  INPUT  "ANOTHER  VALUE?", S* 

530  GOSUB  Check 

540  IF  Check=l  THEN  528 

550  IF  UPC*<S*C1, U)="Y"  THEN  340 

580  INPUT  "ANOTHER  SUBFILE7"^* 

570  GOSUB  Check 

580  IF  Check=l  THEN  560 

580  IF  UPC*<S*Ct, 1])="Y"  THEN  280 

600  GOTO  Menu 

610  Absol!  INPUT  "WHICH  ABSOLUTE  OBSERVATION  NUMBERS", Ob 
620  IF  <Ob>0>  AND  <Ob<«Hobs)  AND  <FRACT<Ob)=0>  THEN  660 
630  PRINT  PAGE;"***  IMPROPER  OBSERVATION  NUMBER  ***" 

640  BEEP 

650  GOTO  Absol 

660  PRINT  PAGE 

670  INPUT  "WHICH  VARIABLE  NUMBER'7" ,  Vn 

680  IF  <Vn>0>  AND  <'Vn<=Nul  AND  < FRACT i  Vn  <=0)  THEN  720 

680  BEEP 

700  PRINT  PAGE;"***  IMPROPER  VARIABLE  NUMBER  ***" 

710  GOTO  670 

720  PRINT  PAGE 

730  DISP  "OLD  VALUE  “;XtOb, Vn>; "  HEW  VALUE  =  “ ; 
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740  INPUT  A 

750  PRINT  USING  760;"0BS  tt:  "  .  Ob,  "  OLD  VALUE:  "  ,  XCOb,  Vn)  ,  "  NEW  VALUE:". A 
760  IMAGE  6A,4D, 12A,5D.5D, 12A,5D.5D 
770  X(Ob, Vn)=A 

7SO  INPUT  "CHANGE  ANOTHER  OBSERVAT I ON^" , ST 

790  GOSUB  Check 

SOO  IF  Check =1  THEN  780 

310  IF  UPC*<S$C1, 1])=“Y"  THEN  Absol 

320  GOTO  Menu 

330  Chgset:  !  PRINT  PAGE 

340  INPUT  "WHICH  SUBFILE  tP",Sn 

350  IF  <Sn>0)  AND  <Sn<=Nsub>  AND  <FRACT'Sn->=0,  THEN  390 
360  BEEP 

370  PRINT  PAGE;  "***SUBFILE  NUMBER  NOT  DEFINED***" 

380  GOTO  340 
890  PRINT  PAGE 

900  PRINT  "WHICH  SET  OF  OBSERVATIONS  IN  SUBFILE  "iNiubt'Sn; 

910  INPUT  Subn 

920  IF  <Subn>0>  AND  <Subn<=Subf <Sn+l )-Subf <Sn  » >  AND  FPACT 1 Subn  1  =0  1  THEN 
930  BEEP 

940  PRINT  PAGE; "***SUBFILE  OBSERVATION  ADDRESS  OUT  OF  PANGE***" 

950  GOTO  900 
960  PRINT  PAGE 

970  PRINT  "OLD  VALUES  FOR  OBSERVATION  ";Subn 

980  PRINT  LIN(2) 

990  Jj=0 

1000  FOR  A=1  TO  Nu 

1010  Jj=Jj+l 

1020  IF  Jj<=5  THEN  GOTO  1050 

1030  Jj  =  l 

1040  PRINT  USING  1060 

1050  PRINT  USING  1070; X< Sub f<Sn>+$ubn, A) 

1060  IMAGE  / 

1070  IMAGE  #,5X,5D.4D 

1080  NEXT  A 

1081  PRINT  LIN<5> 

1090  PRINT  "ENTER  NEW  VALUES  FOR  OBSERVATION  " ; Subn 
1100  PRINT  IIN<2> 

U10  DIM  Aa<20> 

1120  MAT  INPUT  Aa 

1130  FOR  A=  1  TO  Nu 

1140  X<$ubf<Sn)+Subn,  A>=Aa<A.» 

1150  NEXT  A 

1160  INPUT  "ANOTHER  SET  OF  OBSERVAT I ONS , ST 

1170  GOSUB  Check 

1180  IF  Cheek=t  THEN  1160 

1190  IF  UPC*<S*U,  1]>  =  "Y"  THEN  900 

1200  GOTO  Menu 

1210  Add: 1  PRINT  PAGE 

1220  IF  Nsub+ 1 <  =20  THEN  1270 

1230  BEEP 

1240  DISP  “***  MAXIMUM  NUMBER  OF  SUBFILES  EXCEEDED  ***" 

1250  WAIT  3030 
1260  GOTO  Menu 

1270  INPUT  "NUMBER  OF  OBSERVATIONS  IN  SUBFILE''"' ,  No 
1280  IF  Nobs  +  No<  =500  THEN  1320 

1290  PRINT  PAGE;"***  MAXIMUM  HUMBER  OF  OBSERVATIONS  EXCEEDED  ***" 

13O0  BEEP 
1310  GOTO  1270 
1320  PRINT  PAGE 
1330  Nobs=Nobs+No 
1340  Nsub=Nsub+l 

1  350  INPUT  “SUBFILE  NAME0  <5  CHARACTERS  OR  LESS > “ , Nsubf Ulsub * 

1360  Subf  <Nsub+ 1 >=Nobs 
1370  FOR  A= 1  TO  No 
1380  FOR  B= 1  TO  Nu 

1390  DISP  “OBERVAT I  ON  #  ";A;“  VARIABLE  #  ";B;"  IS 

1400  INPUT  X<Nobs-No+A, B> 
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1410  NEXT  B 
1420  NEXT  A 

1430  INPUT  "ANOTHER  SUBFILE9", ST 

144e  GOSUB  Check 

1450  IF  Check=l  THEN  1430 

1460  IF  UPC*<S$U, 13>="V  THEN  1020 

1470  GOTO  Menu 

1480  Addo:  !  INPUT  "DO  YOU  WANT  TO  ADD  AN  OBSERVATION9" , Si 
1490  IF  Nobs+1<=500  THEN  1540 
1500  BEEP 

1510  DISP  "***  MAXIMUM  NUMBER  OF  OBSERVATIONS  EXCEEDED 
1520  WAIT  3900 
1530  GOTO  Menu 

1540  INPUT  "DO  YOU  WANT  SUBFILE  ADDRESSING  C0>  OR  ABSOLUTE  ADDRESSING  S 

1550  IF  <S=0 )  OR  <S=1)  THEN  1590 
1560  BEEP 

1570  PRINT  PAGE;"**-  IMPROPER  RESPONSE  ***“ 

1580  GOTO  1540 
1590  PRINT  PAGE 
1600  IF  S=0  THEN  Suba 

1610  INPUT  "WHICH  ABSOLUTE  OBSERVATION  NUMBER 9 " ,  N 
1620  IF  <N>0)  AND  <H<=Nobi>  AND  <FRACTO1>=0;  THEN  1660 
1630  BEEP 

1640  PRINT  PAGE;"**-  IMPROPER  OBSERVATION  HUMBER  ***" 

1650  GOTO  1610 

1 660  PRINT  PAGE 

1670  FOR  A= 1  TO  Nv 

1680  FOR  B=Nobs  TO  N  STEP  -1 

1S90  X<B+1 , A)=X<B, A) 

1700  NEXT  B 

1710  NEXT  ft 
1720  FOR  A= 1  TO  Nv 

1730  DISP  "INPUT  VARIABLE  #";A;"  FOR  OBSERVATION  I* " ; ft ; 

1740  INPUT  "",Xai,A> 

1750  NEXT  ft 

1760  FOR  A= 1  TO  Nsub+1 

1770  IF  N<*Subf<ft>  THEN  Subf  <A>*$ubf  tAUl 

1780  NEXT  ft 
1790  Nobs=Nobs+l 

1800  INPUT  "ANOTHER  OBSERVATION9" , S* 

1810  GOSUB  Check 

1820  IF  Check3 1  THEN  1800 

1830  IF  UPCi<S*Cl , 1]>=“Y"  THEN  1610 

1840  GOTO  Menu 

1850  Suba:  INPUT  "WHAT  SUBFILE  tt  DO  YOU  WANT  TO  ADD  AN  OBSERVATION  TO9". Sun 
1368  IF  <Sun>0>  AND  <$un<=Njub;»  AND  <FRfiCT<Sun>=0>  THEN  1900 
1870  PRINT  PAGE;"***  IMPROPER  SUBFILE  NUMBER  ♦**" 

1880  BEEP 
1890  GOTO  Suba 
1900  PRINT  PAGE 

1910  INPUT  "WHAT  OBSERVATION  NUMBER  DO  YOU  WANT  TO  ADD9" , Obnuru 
1920  IF  Nobs+t<=500  THEN  1960 
1930  BEEP 

1940  PRINT  PAGE!"***  MAXIMUM  NUMBER  OF  OBSERVATIONS  EXCEEDED  ***" 

1950  GOTO  Menu 
I960  PRINT  PAGE 

1970  IF  <0bnum>0>  AND  •  Obnua.\=$ubf  < Sun*  1  >-Subf  (Sun,* l  >  AND  «  FI  ACT'  Obnuru  =0  •  TH 
N  2010 
1980  BEEP 

1990  PRINT  PAGE;"***  IMPROPER  OBSERVATION  NUMBER  FOR  THIS  SUBFILE  ***" 

2000  GOTO  1910 
2010  PRINT  PAGE 

2020  FOR  A-Nobs  TO  Subf «.Sun)  +Obnuw  STEP  -1 
2030  FOR  B=1  TO  Nu 

2040  X<ft+1 , B)=X<A, B) 

2050  NEXT  B 

2060  NEXT  A 
2070  FOR  ft- l  TO  Nu 
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2030  DISP  "SUBFILE  "  ;  Nsub*<Sun )  ;  “  OBSERVATION  "SObnufti:"  VARIABLE  it " :  A ;  "  IS 

2090  INPUT  X<Obnum+Subf <Sun) , A> 

2100  NEXT  A 
2110  Nobs=Nobs+l 
2120  FOR  A=Sun+l  TO  Nsub+1 
2130  Subf  <  A>=Subt'<A>  +  l 

2140  NEXT  A 

2150  INPUT  "ANOTHER  ADDITION  TO  THIS  SUBFILE7" , ST 

2160  GOSUB  Check 

2170  IF  Check=l  THEN  2150 

2130  IF  UPC*<S*C1, n>="N“  THEN  2200 

2190  GOTO  1910 

2200  INPUT  "ANOTHER  ADDITION  TO  ANOTHER  SUBFILE7". ST 

2210  GOSUB  Check 

2220  IF  Check=l  THEN  2200 

2230  IF  UPCf <STC 1 ,  13)  =  "Y“  THEN  Suba 

2240  GOTO  Menu 

2250  Delete!!  INPUT  "DO  YOU  WANT  TO  DELETE  AN  OBSERVATION7" .Sf 

2260  INPUT  "DO  YOU  WANT  ABSOLUTE  ADDRESSING  'O'  OR  SUBFILE  ADDRESSING  ■  1  7IV: 
2270  IF  <S=0>  OR  <S=1)  THEN  2310 
2280  BEEP 

2290  PRINT  PAGE;"*-*-  IMPROPER  RESPONSE  -  — " 

2300  GOTO  Delete 

2310  PRINT  PAGE 

2320  IF  S=1  THEN  Subdel 

2330  INPUT  "ABSOLUTE  OBSERVATION  NUMBER7" , On 

2340  IF  <On>0>  AND  <On<=Nobs>  AND  (FRACT < On  >=0>  THEN  2330 

2350  PRINT  PAGE;"***  IMPROPER  OBSERVATION  NUMBER  ***" 

2360  BEEP 
2370  GOTO  2330 

2330  PRINT  PAGE;"  OBSERVATION  NUMBER  " ; On 

2390  FOR  A= 1  TO  No 

2400  PRINT  "VARIABLE  tt";A;"  =  ";X<On,A> 

2410  NEXT  A 

2420  INPUT  "DELETE'7"  ,  ST 

2430  GOSUB  Check 

2440  IF  Check=l  THEN  2420 

2450  IF  UPC*<S*Cl,  m  =  “N"  THEN  Abort 

2460  FOR  A=On  TO  Nobs 

2470  FOR  B® 1  TO  No 

2480  X<A,B)=X<A+1,B> 

2490  NEXT  B 

2500  NEXT  A 

2510  FOR  A® 1  TO  Nsub+1 

2520  IF  Subt'<R>  >0n  THEN  Subf < A>»Subf < At- 1 
2530  NEXT  A 
2540  Nobs=Nobs- 1 

2550  Abort:  INFUT  "DELETE  ANOTHER  OBSERVAT ION?" , ST 

2560  GOSUB  Check 

2570  IF  Check  =  1  THEN  Abort 

2580  IF  UPC$<Sf C 1 , 1 ] >=“Y"  THEN  2330 

2590  GOTO  Menu 

2600  Subdel!  INPUT  "SUBFILE  NUMBER  OF  DELETED  OBSERVATION7" , Sn 
2610  IF  <Sn>0>  AND  <Sn<=Nsub>  AND  <FRACT<Snt=0)  THEN  265© 

2620  BEEP 

2630  PRINT  PAGE;"***  IMPPOPER  SUBFILE  NUMBER  *-  +  " 

2640  GOTO  Subdel 
2650  PRINT  PAGE 

2660  INPUT  "OBSERVATION  HUMBER  IN  SUBF I LE7" , On 

2670  IF  i.On>0)  AND  <On<=Subf  tSn+1  >-Subf  tSn;  >  AND  tFRACT*  0n>  =  0^  THEN  2710 
2680  BEEP 

2690  PRINT  PAGE;"***  IMPROPER  OBSERVATION  NUMBER  FOR  THIS  SUBFILE  **-“ 

27O0  GOTO  2660 

2710  PRINT  PAGE;  "SUBFILE:  “ ; NsubT < Sn t ; "  OBSERVATION  it  ";0n 
2720  FOR  A= 1  TO  No 

2730  PRINT  "VARIABLE  tt  ";A;“  IS  “  ;  XtOn+Subf  tSn  » .  A ) 

2740  NEXT  A 
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2750  INPUT  "DELETE?", Si 

2760  GOSUB  Check 

2770  IF  Check=l  THEN  2750 

2780  IF  UPCi<S$Cl, 13)="N"  THEN  Obor 

2790  FOR  A=Subf <Sn)+On  TO  Nobs 

2800  FOR  B=1  TO  No 

2810  X<A,B>=X<A+1,B> 

2820  NEXT  B 

2830  NEXT  A 

2840  FOR  A=1  TO  Nsub+1 

2850  IF  Subf <A>>Subf<Sn)+On  THEN  Subf <A>=Subf <A>-1 

2360  NEXT  A 
2870  Hobs=Nobs-l 

2830  Obor!  INPUT  "ANOTHER  DELETION  FROM  THIS  SUBFILE7" . Si 

2890  GOSUB  Check 

290O  IF  Check=l  THEN  Obor 

2910  IF  UPCi<SiCl, n>="Y"  THEN  2660 

2920  INPUT  "ANOTHER  DELETION  FROM  ANOTHER  SUBFILE7", Si 

2930  GOSUB  Check 

2940  IF  Check=l  THEN  2920 

2950  IF  UPCi<SiCl, 13>="Y"  THEN  Subdel 

2960  GOTO  Menu 

2970  Merge:  INPUT  "FILE  NAME?", Flint 
2930  ON  ERROR  GOTO  Undefined 
2990  ASSIGN  #1  TO  Flint 
300O  READ  Hl;Nv2,Nobs2,Nsub2 
3010  IF  Nobs 2+Nobs<.  =500  THEN  3060 
3020  BEEP 

3030  DISP  "*♦*  MAXIMUM  NUMBER  OF  OBSERVATIONS  EXCEEDED  ***" 
3040  WAIT  2000 
3050  GOTO  Menu 

3060  IF  Nsub2+Nsub<=20  THEN  3110 
3070  BEEP 

3080  DISP  "***  MAXIMUM  NUMBER  OF  SUBFILES  EXCEEDED  ***" 

3090  WAIT  2000 

3100  GOTO  Menu 

3110  IF  Nu2=Nw  THEN  3160 

3120  BEEP 

3130  DISP  "***  VARIABLE  NUMBER  MISMATCH  BETWEEN  MERGED  FILES 

3140  WAIT  2000 

3150  GOTO  Menu 

3160  FOR  A=1  TO  Nsub2 

3170  READ  #1 ; Nsubi<A+Nsub> 

3180  NEXT  A 

3190  READ  tt 1 ; A 

3200  FOR  0=2  TO  Nsub2+1 

3210  READ  #i;Subf<A+Nsub> 

3220  Subf < A+Nsub>=Subf < A+Nsub>+Subf <Nsub+ 1 ) 

3230  NEXT  A 

3240  FOR  A= 1  TO  Nobs2 

3250  FOR  B= 1  TO  20 

3260  READ  # 1 ; X< A+Nobs, B> 

3270  NEXT  B 

3280  NEXT  A 
3290  Nsub=Nsub+Nsub2 
3300  Nobs=Hobs+Nobs2 
3310  OFF  ERROR 
3320  GOTO  Menu 

3330  Undefined:  IF  ERRN=56  THEN  3380 

3340  BEEP 

3350  DISP  ERRMi 

3360  PAUSE 

3370  GOTO  Menu 

3380  BEEP 

3390  PRINT  PAGE;"***  FILENAME  IS  UNDEFINED  ***“ 

3400  GOTO  Merge 
3410  Quit:  ! 

3420  PRINT  PAGE 


39 


APPENDIX  B:  HEWLETT-PACKARD  USERS  MANUALS 

Advanced  Programming  ROM  Manual.  Hewlett  Packard  Co.,  Part  No. 
09845-92065,  Feb.  1980,  108  p. 

BASIC  Language  Interfacing  Concepts.  Hewlett  Packard  Co.,  Part  No. 
09835-90600,  Sept.  1979,  189  p. 

Graphics  ROM  Manual.  Hewlett  Packard  Co.,  Part  No.  09845-91050,  May  1979 
203  p. 

I/O  ROM  Manual.  Hewlett  Packard  Co.,  Part  No.  09845-92060,  Aug.  1980, 

192  p. 

Mass  Storage  ROM  Manual.  Hewlett  Packard  Co.,  Part  No.  09845-92070,  Feb. 
1980,  132  p. 

System  45  Operating  and  Programming  Manual.  Hewlett  Packard  Co.,  Part  No 
09845-92000,  Feb.  1980,  302  p. 
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